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1.0 BACKGROUND 


As detailed in the original Scope of Work for this research 
effort, the principal investigator was to improve the steady-state 
power balance model (PBM) of the Space Shuttle Main Engine (SSME) 
in a three phase effort. A summary of the tasks in each phase is 
given below. 

Phase 1: Construct software to facilitate SSME performance 
prediction and test data validation. 

Phase 2 : Review computational logic within the current 
version of the SSME PBM and implement programming structure. 

Phase 3 : Develop programming logic to improve the physical 
consistency of routing routines within the SSME PBM. 

After discussions with John Butas of the NASA/MSFC Propulsion 
Laboratory in January of 1992, the contract Scope of Work for Phase 
3 was modified to place greater emphasis on refinement of 
computational tools initially developed in Phases 1 and 2 , and to 
support evaluation of recently acquired Technology Test Bed (TTB) 
data. The primary effort during Phase 3 was continued development 
of software intended to support integration of TTB test data with 
SSME performance predictions from the steady-state power balance 
model. In addition, variational analyses of both TTB test data and 
PBM predictions were performed. Results of these analyses were 
compared to evaluate the computational integrity of the SSME 
steady-state power balance model. 

A considerable portion of the contract effort was dedicated to 
development and testing of a formal strategy for reconciling 
uncertain test data with physically limited computational 



prediction. This emphasis was motivated by the availability of an 
extensive and highly organized SSME performance data base from the 
Engine 3001 test program, and by serious inconsistencies in power 
balance model predictions. 

A review of SSME steady-state power balance model function is 
provided in Section 2 of this report. Specific weaknesses in the 
logical structure of the current PBM version are described with 
emphasis given to the main routing subroutines BAL and DATRED. 
Selected results from a variational analysis of PBM predictions are 

compared to TTB variational study results to assess PBM predictive 
capability. 

The motivation for systematic integration of uncertain test 
data with computational predictions based on limited physical 
models is provided in Section 3 . The theoretical foundation for 
the reconciliation strategy developed in this effort is presented, 
and results of a reconciliation analysis of the SSME high pressure 
fuel side turbopump subsystem are examined. Specific 
recommendations are presented in Section 4 . 



2.0 SSME PBM LOGIC ASSESSMENT 


The Space Shuttle Main Engine power balance model is a FORTRAN 
based software package developed by the Rocketdyne Division of 
Rockwell International. it is used to predict operating 
characteristics and performance of the SSME under steady-state 
conditions. Approximately 800 SSME temperatures, pressures, flow 
rates, shaft speeds, and other hardware performance parameters are 
calculated by the power balance model. 

The current version of PBM has a number of analysis options. 
The standard power balance analysis option determines fluid and 
flow properties throughout the entire engine system assuming 
nominal hardware performance characteristics. In addition, there 
is a data reduction analysis option which uses actual test data to 
define hardware operating characteristics such as efficiencies and 
flow multipliers of a specific SSME. The base balance option is 
used to further define engine hardware characteristics by matching 
nine parameters to data reduction output. 

Although conceptually a powerful prediction tool, PBM exhibits 
a number of significant shortcomings. Documentation of the 
computational, physical, and functional operation of PBM has not 
been rigorously maintained and is inadequate. Moreover, recent 
tests have demonstrated that PBM predictions fail to satisfy 
fundamental energy balance relations within all engine subsystems 
[1]. As a result, confidence in PBM predictions has been degraded 
and software utility diminished. 

In order to assess the logical integrity of PBM, five sources 



of information were utilized: 


1) detailed flowcharts of the main analysis routing routines BAL 
and DATRED 

2) iteration loop diagrams for routines BAL and DATRED 

3) detailed translation of routines BAL and DATRED 

4) direct source code inspection 

5) analysis of variations (ANOVA) comparisons of PBM predictions 
and TTB data. 

An overview of the main routing logic was obtained from sources 1 
and 2. The physical and logical consistency of individual lines of 
code was examined using sources 3 and 4. The integrity of the PBM 
prediction process was evaluated using comparisons from source 5. 

Detailed, computer-generated flowcharts of PBM subroutines BAL 
and DATRED were obtained from NASA/MSFC/EP52 and examined for 
logical structure. Copies of these flowcharts were previously 
presented in the Phase 2 final report. A high level of vertical 
connectivity (logic feedback) is obvious upon examination of the 
complex BAL flowchart. Subroutine logic is highly integrated and 
structured segmentation cannot be achieved without fundamental and 
costly code modifications beyond the scope of this effort, 
^^^routine DATRED has a more sequential logic process, however, 
structured segmentation of the code was not attempted for reasons 
described below. 

Iteration loop logic diagrams for subroutines BAL and DATRED 
were also presented in the Phase 2 final report. Subroutine BAL 
contains four multivariate iteration loops for solution of 



simultaneous nonlinear relations, and twenty-five univariate 
iteration loops for solving individual nonlinear relations, all of 
the Newton-Raphson type. Five deep nesting of iteration sequences 
is found in BAL with high level multivariate iteration loops 
traversing virtually the entire routine. Intersection of major 
iteration sequences inhibits structuring of code logic. Iteration 
nesting and crossover are not as severe in subroutine DATRED, 
however, both BAL and DATRED use a segmented solution strategy on 
restricted subsets of the fluid and flow governing equations. 
Values of the subset solution variables are iteratively matched 
with both nested and sequential subset solutions. This type of 
segmented solution approach with matching is generally less 
than robust global strategies for solution of nonlinear 
equations [see, e.g. 2]. 

In order to facilitate interpretation of PBM logic, a software 
^' rans l a ^- or package was constructed in the C programming language. 
The translator substitutes variable definitions found in PBM 
documentation in place of the variable names in PBM. The result is 
a readable document describing PBM function in detail. Translations 
of subroutines BAL and DATRED were included in the Phase 2 final 
report. These translations were used to study software logic line 
by line. 

A detailed examination of BAL and DATRED logic indicates that 
many computations are empirical and/or heuristic. This conclusion 
is based on comparison of the actual number of SSME flow network 
controllers with the number of independent variables used by PBM to 
predict a variety of engine operating conditions. Since the SSME is 



a feedback dominated flow network, each control setting can be 
expected to affect operating characteristics throughout the engine. 
However, many PBM computations are based on reduced dependencies. 
This is especially evident in the data reduction routine DATRED 
where certain densities, temperatures, flow rates, pressures, and 
hardware characteristics are specified by relations depending on 
thrust level (or commanded chamber pressure) alone. 

In order to assess the fundamental dependencies of FBM 
computations, a variational study of power balance predictions was 
performed. PBM analyses of engine number 3001 were performed by 
John Butas of the MSFC Propulsion Laboratory. Analysis independent 
parameters were set at values corresponding to control parameter 
se htings for each of the 16 test profiles employed in the TTB 
Engine 3001 test program. Control parameter definitions and 
settings for each of the 16 TTB program tests are shown in Appendix 
B, Table 1. 

The Engine 3 001 test series was based on a Taguchi type design 
of experiments [ 3 ] . The matrix of control settings displayed in 
Table 1 was selected based on a fractional factorial test plan to 
facilitate data utilization. A variational analysis of TTB 
recorded engine operating conditions at test matrix control 
settings was also performed. Results of the TTB data variational 
s hudy were then compared to the PBM analysis of variations. 

parameter comparisons are displayed in Appendix A, Figures 
^ through 4 . The computed contribution of each individual control 
parameter to the variation of the performance variable listed in 
the title is displayed. The abscissa designations OSP, HSP, ORP, 



and HRP correspond to LOX NPSP, FUEL NPSP f LOX REPRS , and FUEL 
REPRS respectively in the test matrix. The category COMB that 
appears in each figure represents contributions from control 
parameter combinations that cannot be individually allocated 
because the test program was not designed as a full factorial set 
of experiments. 

It is obvious from Figure 1 that, within the TTB test range, 
low pressure fuel pump (LPFP) speed variation was only weakly 
affected by F7 orifice size. pbm predicted F7 contributions to 
LPFP speed variation were significantly greater than test results 
indicated. Similarly, as shown in Figure 2, TTB data indicated a 
significant F7 orifice size contribution (17%) to high pressure 
fuel pump (HPFP) discharge temperature variation that was largely 
absent from PBM predictions. Comparisons such as these are 
indicative of potential PBM weakness in assigning component level 
contributions to performance and in predicting operational 
contributions of hardware modification. 

f ^^snces (approximately 31%) between observed and 
predicted power level (% RPL) and mixture ratio (M/R) contributions 
to coolant control valve (CCV) flow rate are indicated in Figure 3. 
Significant disparity between predicted and observed controller 
effects was not isolated to the parameters displayed in Figures 1 
through 3 . In general there was good agreement between predicted 
and observed controller contributions to pressure variation such as 
is shown in Figure 4 for the high pressure fuel pump (HPFP) outlet 

pressure. Only isolated cases of significant pressure variational 
differences were observed. 


More common were large differences 



between predicted and observed contributions to temperature, flow 
rate, and hardware performance characteristics. 

The variational analysis comparisons described above suggest 
that PBM does not adequately model the variation of important SSME 
performance parameters. This is not particularly surprising since, 
in many places throughout the code, physical dependencies have been 
replaced by "hardcoded" empirical relations that lack adequate 
documentation to assess application and range validity. These 
comparisons also suggest that the power balance model would be 
inadequate as a design or anomaly resolution tool. Integrity of 
PBM predictions can be expected only in nominal engine operating 
ranges over which code empirical relations were established. 

Because the power balance model is a highly connected software 
package with significant iteration looping, it is difficult to 
access the overall impact of an individual code modification 
without significant computational testing. Simple code corrections 
to achieve improvements in isolated parameter prediction can have 
a far reaching and detrimental affect. Therefore, code maintenance 
and modification time will be substantially greater than for a 
highly structured, modular, and well documented performance model. 

One of the major functions of the power balance model is to 
integrate test data with theoretical predictions. The weaknesses 
of the existing data integration procedure will be discussed in the 

following section, and a new integration strategy will be 
introduced. 



3.0 RECONCILIATION MODEL 


One of the features of the steady-state power balance model is 
its ability to integrate test data into the performance prediction 
process. This is accomplished within the data reduction analysis 
option. In the data reduction process, test information is 
incorporated literally into predictions, and hardware operating 
parameters are adjusted to values consistent with this presumed 
pristine data. 

Unfortunately, experimental data associated with a complex 
flow system such as the SSME is fraught with uncertainty. 
Maintaining operation and calibration of sensing and signal 
conditioning instrumentation is difficult in the severe SSME 
operating environment. In addition, point measurements in such a 
complex flow environment often include the effects of highly 
localized and/or secondary flow phenomena that are not 
characteristic of system average conditions. Literal incorporation 
of inaccurate test data can lead to nonphysical predictions of 
engine operation and erroneous assumptions concerning hardware 
• Since all test data has associated uncertainty, the 
pristine data assumption is inappropriate in a test information 
integration strategy. 

Performance prediction models based on fundamental flow 
physics are also limited by theoretical approximations required to 
achieve tractable solution. For example, PBM computations assume 
steady-state operation throughout the engine, and provide estimates 
of average flow conditions using a cross-stream uniform, one- 



dimensional flow approximation. In addition, thermodynamic 
property data for hydrogen and oxygen in SSME operating ranges has 
accuracy limitations. These type assumptions and limitations 
necessarily restrict the accuracy of theoretical model predictions 
and present an additional source of uncertainty for data 
integration strategies. 

The above observations suggest that any systematic procedure 
for integrating experimental data and theoretical predictions 
should recognize both data uncertainty and model limitations. The 
objective of the reconciliation development effort undertaken as 
part of this study was to construct a logical strategy for 
integrating uncertain test data with limited theoretical 
predictions in order to determine most plausible SSME operating 
conditions. 

A heuristic yet logical procedure for achieving systematic 
data integration was presented in the Phase 1 final report of this 
study. A refinement of this reconciliation procedure has been 
developed in Phase 3 . The basis of the new method rests on the 
principle that the mean of experimental observations reflects most 
probable, but not absolute, engine operating conditions. If 
measured engine operating properties are assumed to be independent, 
normally distributed, random variables, then the most probable 
engine state will maximize the property joint probability density 
function (pdf) subject to constraints imposed by physical laws. A 
mathematical expression for this state pdf is given below: 



F( X 1 ... X k ) 



a, ... a k ( 2?r) k/2 


where 



X, - adjusted value of property i 
<7; ~ standard deviation of property i 

Mf - mean of property i 

d. - deviation of adjusted property i from its mean 
(measured value) ( d, = X { - n . ) 

F - joint probability density function of state 

k - number of properties 

(X 1 . . . X k ) - state of system 


Properties in the relation above include measured flow rates, 
temperatures, and pressures throughout the engine system. 

The state pdf is a maximum when the expression in brackets in 
the exponent of e is minimized. In the absence of physical 
constraints this minimum would occur when all the d ( are zero, or 
when the value of each property i is at its mean. Therefore, 
experimental property measurements are assumed to correspond to the 
property means in equation 1, and the d { are adjustments from 
measurement values required to adequately satisfy physical 
constraints, including mass and energy conservation requirements as 
well as second law limits. 

A robust data reconciliation strategy must also incorporate 
measurement system uncertainty limits in addition to physical 
constraints. The problem of determining most plausible engine 
operating conditions can thus be reduced to a mathematical 
programming problem of the form: 



maximize Z = F(d 1 ... d k ) by selection of (d, ... d k ) 

assuming constant (a 1 ... <j k ) 

subject to 

physical constraints for each engine subsystem j 
| mass flow imbalance ■ | < L flow _j 

| energy imbalance . | < 1^, 

| entropy production . | >0 


uncertainty limits for measurements at each node i 


j mass flow adjustment . 

1 

U nM 


| pressure adjustment . 

1 

U p-i 


| temperature adjustment 1 

1 

U T-i 

(2) 
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mass flow adjustment at node i 

- pressure adjustment at node i 
temperature adjustment at node i 

- upper limit of flow imbalance for flow circuit j 

- upper limit of energy imbalance for control 
volume j 

- flow uncertainty limit at node i 

- pressure uncertainty limit at node i 
temperature uncertainty limit at node i 

” mass flow rate across inlet or outlet surface 

- specific enthalpy at inlet or outlet surface 
specific entropy at inlet or outlet surface 
mass density at inlet or outlet surface 

~ cross-sectional area of inlet or outlet surface 

- rate of energy transfer as heat into control 
volume 

- rate of energy transfer as work from control 
volume 

- temperature at which energy as heat enters volume 


The mathematical programming problem stated in formulation 2 
above is highly nonlinear. Without loss of generality, the 
objective function Z = F can be replaced by the exponent of e in 
equation 1. If, in addition, the imbalance relations are 
approximated as first order truncated Taylor series expansions in 
the nodal adjustment values d, the mathematical programming problem 
reduces to the following: 


k 


minimize Z 


E 

i=l 



k = number of 
measurements 


subject to 

linearized forms of the physical constraints for each 
engine subsystem j 

| LlmbMf a )j | < L floy ,j 

| LlmbE ( a ), | < ^Vwgy-j 

| Limbs ( a )j | > o 


measurement uncertainty limits for each node i 

(n = number of nodes = number of measurements/ 3 = k/3) 
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The objective function Z in formulation 3 above is quadratic 
in the measurement adjustments d, and the constraints have been 
linearized in terms of the adjustment variables d. This is the 
form of a classical quadratic programming problem for which a 
variety of robust solution strategies exist. The solution of this 



problem minimizes the property adjustments required to satisfy 
physical constraints within measurement system uncertainty bounds, 
and in a logical sense provides most plausible engine operating 
conditions subject to restrictions inherent in the linearization of 
the physical constraints. 

The reconciler developed as part of this effort constructs the 
quadratic programming problem defined in formulation 3 above and 
implements the complementary pivot method algorithm [4] to obtain 
the QP problem solution. a hierarchy diagram describing the 
organization of routines in the reconciler is presented in Appendix 
A, Figure 5. Documentation describing the function of reconciler 
routines is given in Appendix Cl, and a source code listing is 
presented in Appendix C2. 

In order to perform a reconciliation analysis, four types of 
input data are required. Thermodynamic property data in operating 
ranges of interest are necessary. Specific enthalpy, specific 
entropy, and density as functions of pressure and temperature are 
required. For SSME analyses, hydrogen, oxygen, and water property 
information was provided and integrated into the reconciler logic 
by John Butas of NASA/MS FC/EP5 2 . In addition, experimental data 
(or computational simulation results) are required to provide a 
baseline for adjustment calculations. The TTB Engine 3001 test 
program has provided extensive high quality experimental data for 
reconciliation analyses. PBM predictions have provided a 
simulation baseline for initial reconciler testing. The third type 
of input required for reconciliation analysis is uncertainty 
estimates quantifying model limitations (physical constraint bounds 



L in formulation 3) as well as test data confidence bands 
(uncertainty bounds U in formulation 3). Finally, system 
definition information must be constructed to specify engine 
configuration and to properly associate nodal property data. a 
detailed description of input data requirements is provided in 
Appendix Cl documentation and a listing of the reconciler input 
data file format is provided in Appendix B, Table 2. 

Reconciler performance has been verified on a number of test 
problems. Recently, reconciler logic was tested using results of 
a PBM simulation of the HPFTP system at 109% RPL to provide 
baseline measurements. A schematic of the HPFTP system with 
analysis nodes identified is displayed in Appendix A, Figure 6. 
The analysis configuration was composed of 14 nodes, 5 mass flow 
circuits, and 4 energy volumes. Mass circuit and energy volume 
definition nodes for this analysis are specified in Appendix B, 
Table 3. The energy volumes include the fuel preburner, high 
pressure fuel turbopump, fuel side turn around duct, and fuel side 
hot gas manifold. Coarse measurement system uncertainty estimates 
were utilized in the HPFTP test case analysis because more precise 
uncertainty information was unavailable. These estimates are 
provided in Appendix B, Table 4. 

Mass, energy, entropy, and availability imbalances both before 
and after reconciliation are displayed in Appendix B, Tables 5 and 
6 respectively. a significant reduction in subsystem energy 
imbalances after reconciliation is the most striking result 
observed in comparing Tables 5 and 6. A 99% energy imbalance 
reduction in the fuel preburner and turbopump subsystems 


was 



obtained during the reconciliation process while mass balance and 
entropy production requirements were maintained. System properties 
both before and after reconciliation are presented in Appendix B, 
Table 7. The adjustments required to achieve reconciliation 
(solution to the quadratic programming problem outlined in 
formulation 3) are also presented in Table 7. Significant 
reductions in PBM predicted hot gas temperatures throughout the 
system are suggested. These temperature reductions remain within 
specified measurement uncertainty bounds, yet provide substantial 
improvement in the initially large energy imbalances. 

Heuristic data integration procedures do not provide the level 
of confidence in prediction that is required in a long term engine 
development program. Efficient and reliable use of experimental 
observation to improve performance prediction and safety requires 
a systematic test data integration strategy. The reconciliation 
strategy outlined above is a logical procedure for improving rocket 
engine performance prediction. The mathematical foundation is well 
defined and computations are physically sound within approximation 
limits. In addition, the base procedure is completely general, 
with material and system configuration provided by modular data 
file inputs. Initial test results have been quite successful and 
strongly support continued development and use of this mathematical 
programming approach for large system test data reconciliation. 
This technique can also be utilized to evaluate test data integrity 
and isolate measurement system problems. 



4.0 RECOMMENDATIONS 


A list of recommendations based on results of this research 
effort is presented below. 

1. Local modifications to the power balance model should be 
thoroughly investigated before implementation due to the high level 
of logic connectivity. If PBM is to be used extensively as a 
prediction tool, a catalog of parameter influence coefficients 
should be developed to efficiently assess the system wide impact of 
specific code changes. 

2. Without extensive documentation describing imbedded 
empiricisms within PBM logic, the power balance model should be 
considered a high order gains model containing the experience base 
archive. PBM should not be considered a cornerstone theoretical 
prediction tool without major modification. 

3 . Development of mathematical programming approaches to test 
data reconciliation should continue in order to provide a 
consistent and logical basis for improving performance prediction, 
a platform for logically resolving data inconsistencies, and a 
means of assessing data and measurement system integrity. 

4. A fundamentally sound theoretical model of engine 
performance should be developed. 

5. Uncertainty analysis should be incorporated in any rocket 
engine performance evaluation and prediction program. 

6. An integrated rocket engine performance prediction 
platform should be developed that modularizes fundamental 
theoretical computations and provides a standardized interface for 



efficient parametric integration of test data. 

7. Frictional resistance relations should be added to the 
quadratic reconciler in order to provide more consistent pressure 
loss relations. 
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FIGURES 


Figure 1 . LPFP Speed 
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Figure 2. HPFP Disch Temp 
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Figure 3. CCV Flow Rate 
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Figure 4. HPFP Disch Pressure 
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Figure 5. Quadratic Reconciler Hierarchy Diagram 
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Figure 6. HPFTP system with reconciliation analysis nodes indicated 
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TABLES 


Table 1. TTB Engine 3001 test program matrix of control variable settings 
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Table 2. Reconciler input data file formats 


FILEIO.DAT - designates I/O data filenames 

'input property data filename' 

'input test data filename' 

' input uncertainty estimates filename ' 
'input volume definition data filename' 
'output (standard format) filename' 

'output (test data input format) filename' 


TDAT = 'input test data filename' 


NDESC NTTB 

DESC(l) DESC (NDESC) 

TTB(l) TTB(NTTB) 


UDAT = 'input uncertainty estimates filename 1 


IPRPD ITTBD IROWA JCOLA ITORDX ITORDY DPF DTF 

UP(1) UP (NTNOD) 

UT(1) UT (NTNOD) 

UW(1) UW (NTNOD) 

UWMFC(I) UWMFC (NTMFC) 

UEVOL(l) UEVOL(NTVOL) 

USVOL(l) USVOL (NTVOL) 


(continued next page) 


/' 

Table 2. Reconciler input data file formats 


VDAT - 'input volume definition data filename' 


IENV IPCTTH 

IA(1) 

IP(1) 

IT(1) 

IW(1) 

MAT ( 1) 

MI0(1) 

NIO(l) 

MODIR (1,1) 

IMFCN (1,1) 


MAXSTG NHG NTMFC NTNOD NTVOL 
• • • • XA (NTNOD) 

IP (NTNOD) 

IT (NTNOD) 

IW( NTNOD) 

MAT (NTNOD) 

MIO (NTMFC) 

NIO (NTVOL) 

MODIR(l,MIO(l) ) 

IMFCN (1, NTMFC) ) 


MODIR (NTMFC, 1) MODIR ( NTMFC , MIO ( NTMFC) ) 

Tr,nxT 3 ^ F ^ (N ™ FC,1) IMFCN (NTMFC, MIO (NTMFC) ) 

IODIR (1,1) IODIR ( 1 , NIO ( 1 ) ) 

. IVOLN (1,1) IVOLN ( 1 , NTVOL) ) 


IODIR(NTVOL, 1) IODIR (NTVOL, NIO (NTVOL) ) 

Mon H r^T? LN(NTV ° L,1) IVOLN (NTVOL, NIO (NTVOL) ) 

NODHG ( 1 ) NODHG ( NHG ) ' ' 

NH2HG ( 1) NH2HG (NHG) 

N02HG ( 1) N02HG (NHG) 

ICEFF(l) ICEFF (NHG) 

IH2HG (1,1) IH2HG ( 1 , NH2HG ( 1) ) 

I02HG (1,1) I02HG ( 1 , N02HG ( 1 ) ) 


IH2HG (NHG, 1) IH2HG (NHG, NH2HG (NHG) ) * 

102 HG (NHG, 1) I02HG (NHG , N02HG (NHG) ) 



Table 3 . HPFTP analysis number 1 circuit definitions 


Energy volume Definitions 
Volume # Boundary Nodes 

1 12 3 4 

2 4 5 6 7 

2 6 7 8 9 

4 9 10 11 12 


Mass Circuit Definitions 


Circuit # 
1 
2 

3 

4 

5 


Boundary Nodes 

12 3 4 

5 6 14 

4 7 14 

6 8 13 

7 9 13 



Table 4. Uncertainty estimates for HPFTP analysis number l 


Nodal Property Uncertainty 


Node # 

Pressure 

(psi) 

1 

0.772 

2 

0.605 

3 

0.659 

4 

200.0 

5 

0.0272 

6 

200.0 

7 

200.0 

8 

200.0 

9 

200.0 

10 

200.0 

11 

0.359 

12 

200.0 

13 

200.0 

14 

200.0 


Temperature 

Flow Rate 

(degF) 

(lb/sec) 

0.0208 

0.00796 

0.1 

0.00831 

0.00972 

0.0001 

100.0 

1.6 

0.00424 

0.0162 

20.0 

1.7 

100.0 

1.6 

20.0 

1.6 

100.0 

1.7 

100.0 

1.7 

0.0436 

0.00175 

200.0 

0.00175 

20.0 

0.14 

20.0 

0.3 


Hass Imbalance Limits 


Circuit # 

Flow Imbalance Limit 
(lb/sec) 

1 

0.001 

2 

0.001 

3 

0.001 

4 

0.001 

5 

0.001 


Volume # 

Volume Imbalance Limits 
Energy imbalance Limit 

Entropy imbalance 


(Btu/sec) 

(Btu/degR-sec) 

1 

71.0 

100.0 

2 

43.0 

100.0 

3 

2.6 

100.0 

4 

18.0 

100.0 


Table 5 . High pressure fuel turbopump system imbalances at 
109% RPL prior to reconciliation 


HP_FUEL_S I DE_ANALYSIS_1 
THRUST « 109.0% RPL 


SUBSYSTEM DW(LB/S) DE(BTU/S) DE (HP) DS(BTU/R-S) DAV(BTU/S) 


FPB_SUBSYSTEM 
HPFTP_SUBSYSTEM 
TAD_SUBSYSTEM 
HGM SUBSYSTEM 


.014 -7111.634 
.000 -4277.324 
.100 -2563.647 
.000 183.112 


10061.800 -664 
-6051.713 -174 
-3627.140 -12 

259.074 -5 


.165 

344895.60 

.517 

88216.79 

.175 

3889.41 

.406 

3048.38 


Table 6. High pressure fuel turbopump system imbalances at 
109% RPL after reconciliation 


HP_FUEL_S I DE_ANALYS IS_1 
THRUST = 109.0% RPL 


SUBSYSTEM 


DW(LB/S) DE (BTU/S) 


DE (HP) DS ( BTU/R-S ) DAV( BTU/S) 


FPB_SUBS Y STEM 
HPFTP_SUBSYSTEM 
TAD_SUBSYSTEM 
HGM SUBSYSTEM 


.002 

-59.775 

-84.571 

-660.372 

.003 

-17.536 

-24.811 

-165.596 

.004 

-58.564 

-82.858 

-6.975 

.024 

-227.164 

-321.400 

-.467 


349937.30 

87748.05 

3638.13 

20.30 



Table 7 


HPFTP System Reconciliation at 109% RPL 


HPFTP ANALYSIS NO 1 


ORIGINAL NODE DATA 


>DE 

PRESSURE 

1 

7723.00 

2 

6049.00 

3 

6585.00 

4 

5498.00 

5 

271.70 

6 

6739.00 

7 

3718.00 

8 

6739.00 

9 

3691.00 

10 

3585.00 

11 

3587.00 

12 

3574.00 

13 

6739.00 

14 

6739.00 


TEMPERATURE 

flowrate 

208.30 

79.63 

278.40 

83.07 

97.16 

1.01 

1929.00 

163.70 

42.40 

162.30 

96.71 

159.20 

1737.00 

166.80 

96.71 

157.80 

1720.00 

168.10 

1715.00 

168.10 

436.30 

17.51 

460.90 

17.51 

96.71 

1.38 

96.71 

3.05 


RECONCILED NODE DATA 


NODE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


PRESSURE 

7723.00 

6049.00 

6585.00 
5511.55 

271.70 

6719.12 

3704.84 

6723.15 

3657.60 

3569.42 

3587.00 
3543.89 
6722.71 
6722.71 


TEMPERATURE 

208.30 
278.40 

97.16 

1907.02 
42.40 
95.07 

1704.95 

93.87 

1682.59 

1687.03 

436.30 
419.48 

96.53 

96.53 


FLOWRATE 

79.63 

83.07 

1.01 

163.71 

162.30 

159.25 

166.76 

157.87 

168.14 

168.12 

17.51 

17.51 

1.38 

3.05 


(continued next page) 


Table 7. HPFTP System 


at 109% RPL 



.00 

.00 

.00 

-21.98 

.00 

-1.64 

-32.05 

-2.84 

-37.41 

-27.97 

.00 

-41.42 

-.18 

-.18 


.00 
.00 
. 00 
.01 
.00 
.05 
-.04 
.07 
.04 
.02 
.00 
.00 
.00 
.00 


.00 

.00 

.00 

-1.14 

.00 

-1.69 

-1.84 

-2.94 

-2.18 

-1.63 

.00 

-8.99 

-.18 

-.18 




APPENDIX Cl 


RECONCILER DOCUMENTATION 



Routine: MAIN 


(main program 


for reconciliation model) 


Routine Function: 

M . „ aVMna routine for reconciliation model. Opens and 
Mam cal11 ™? ^ ssoci ated with thermodynamic 

uJooJrtHata lon«fl volSme dlfinition, test data, and 
uncertainty^estimates . Initializaes parameters for 
reconciliation model • 


Common Blocks: 

Blank - 

TDAT - 
UDAT - 
VDAT ~ 
H2PRP - 

02PRP - 

H20PRP- 

TABLE - 

STD 


quadratic programming algorithm matrices 
and parameters 
test data 

uncertainty estimates . . ■ 

control volume definition informatio 
hvdroqen table pressures and temperatures, ... 

and lUSciated Snthalpies, entropies, and densities 
oxvaen table pressures and temperatures, . . 

and 9 associated enthalpies, entropies, and densities 
water table pressures and temperatures, . . 

andassociated enthalpies, entropies, and densities 
number of distinct pressures and temperatures 
each hydrogen, oxygen, and water table absolute 
reference state enthalpies entropies, and absolute 
entropies for hydrogen, oxygen, and water 


Input Variable Definitions: 


Input File 
PDAT 

TDAT 

UDAT 

VDAT 

ODAT 

ODATV 


' FI1 f I alpha ' variable'that 1 identifies’ name of property 

- alpha 1 variabie S that identifies name of test data 

input file - 

- alpha variable that identifies name of 

uncertainty data input fil® , rae 

- alpha variable that identifies name of volume 

defintion input file 

- alpha variable that identifies name of 
main reconciler output fil® 

- alpha variable that identifies name of 
reconciler output file that is m same format as 
variable TDAT input file 


Input File PDAT 

PTITLE 
NH2P (I) 
NH2T ( I ) 


(fil^nimc^variable^identified in FILEIO.DAT, 

- property data table subsection title . T 

- number of pressures in H2 property data table 

- number of temperatures in H2 property data 
table I 


Routine MAIN 


(page 2) 


H2P1 ( I ) 

H2T1 ( I ) 
H2H1(I,J) " 

H2S 1 ( I / J ) _ 

H2D1 ( I / J) " 

H2P2 (I) 

H2T2 (I) 

H2H2 (I r J) - 

H2S2(I,J) " 

H2D2 ( I , J) “ 

H2P3 ( I ) 

H2T3 ( I) 

H2H3 ( I , J) - 

H2S3 ( I , J) - 

H2D3(I,J) ‘ 

H2P4 (I) 
H2T4 (I) 
H2H4 (I, J) 

H2S4 (I, J) 

H2D4 (I , J)' 

N02P ( I ) 
N02T ( I ) 

02P1 ( I ) 
02T1 (I) 
02H1 ( I , J) 

02S1(I, J) 

02D1(I, J) 

02P2 (I) 
02T2 (I) 
02H2 (I, J) 

02S2(I,J) 

02D2 ( I / J) 


and 


I and 


and 


I and 


pressure I in H2 data table 1 
temperature I in H2 data table 1 and 

H2 enthalpy associated with pres 

temperature J in H2 data J fissure I 
H2 entropy associated with P 
temperature J in H2 data ba ^! s i ure 
H2 density associated with pre 
temperature J in H2 data ta 
pressure I in H2 data tabl 

temperature I in H2 , d p fJ^ Assure I and 
H2 enthalpy associated with pr< ess 

temperature J in H2 data bab *! s J ure I 
H2 entropy associated with press 
temperature J in H2 data nressure 
H2 density associated with press 
temperature J in f data f b 6 
pressure I in H2 data table 3 

temperature I in H 2 d ata Jf^essure I and 

• H2- enthalpy associated with pressure 

temperature J in H2 data Pressure I and 
■ H2 entropy associated with p 
temperature J in H2 data table 

- H2 density associated with pressure 

. sesttiS s 

: srssss JiSsrr " 1 and 

- i -a 

temperature J in H2 data table 4 

- H2 density associated with pressure 

temperature J in H2 data table 4 j 

- number of pressures in 02 property aa 

- number of temperatures in 02 property aa 

table I . , , , 

- pressure I in 02 data table 1 

I S-SKSS »» “atS^ig pressure X and 

- SHn^r-siSia^tn essure X and 

. srs^^-KsrH 1 and 

temperature J in 02 daba Jr ab * e 1 

- pressure I in 02 data table 2 

I o^enthalpy associated with^ressure I and 

- rrntroprassoSiated^lth pressure I and 

- ia^H S H 1 and 

temperature J in 02 data table 


Routine MAIN 


(page 3) 


02P3 ( I ) 

02T3 ( I ) 

02H3 ( I / J) " 

02S3(I,J) “ 

02D3(I,J) “ 

NH20P (I) - 

NH20T (I) 

H20P1 ( I ) 
H20T1 (I) “ 

H20H1 (I » J) ” 

H20S1 ( I r J) “ 

H20D1 ( I , J) “ 


Input File TDAT 

NDESC 

NTTB 

DESC(I) 

TTB(I) 

Input File VDAT 

IENV 

IPCTTH 

MAXSTG 

NHG 

NTMFC 

NTNOD 

NTVOL 

IA(I) 

IP(I) 

IT (I) 

IW(I) 

MAT (I) 

MIO(I) 


pressure I in 02 data table 3 ^ 

rintb^py associated^witlvpressure I and 

ol^entropy^ssociate^with ^pressure I and 

5 rrenfity e ass-iated a with a pressure I and 

number 3 of^ressures in f O property data table 1 
number of temperatures in H20 prope y 

pressure I in H20 data L 

temperature I in H2 a pressure I and 

H20 enthalpy associated with PJ sssur 

temperature J in H2 °^ t h pressure I and 
H20 entropy assoc ^ e * l 

temperature J m Assure I and 

H20 density assoc ^ e * ? 

temperature J in H20 data table 

(test data input ^fdentif ied in FILEI0.DAT) 
(file name vanabl test data descriptions 

. number of alpha dala table 

* number of data entries i 

- alpha’ variable test data description I 

- test data entry I 

(volume definition i n FILEIO.DAT) 

. VST- 

- TTB array address of % US ed in SQP iteration 

- number of stages to be used m =>u 

: = “ zgrzsr 

analysis engine system analysis 

- number of nodes xn^ 9 . svstem analysis 

- number of volumes m 9 t • inq the value of 

- position in TTB array containing tne 

the cross-sectional area f aining the value of 

- position in TTB array containing tn 

- pos it ion* it^TTB S?rV containing the value of 
I “LtloHn^B frrafcontaining the value of 

- material IdentlSlnfn^r at node I 

- JmL>er'of i/0'i associated 5 with mass flow 
circuit I 



Routine MAIN 


(page 4) 


NIO(I) 
MODIR(I, J) - 

IMFC ( I , J ) - 

IODIR(I, J) - 

IVOLN ( I , J) - 
NODHG(I) “ 
NH2HG(I) - 
N02HG (I) - 

ICEFF(I) - 

IH2HG ( I , J) “ 
I02HG(I , J) - 


Input File UDAT 

IPRPD 

ITTBD 

IROWA 

JCOLA 

ITORDX 

ITORDY 

DPF 


DTF 


UP(I) 

UT (I) 
UW(I) 
UWMFC(I) 

UEVOL(I) 

USVOL(I) 


, j: -t /<-) i associated with volume I 

number of I/O s associ circuit I 

flow direction of I/O J J- n “Jr* 

1 = inlet flow, -1 = outlet flow 

node number of I/O J in mass flow circuit 

flow direction of I/O J in volume I 

1 = inlet flow, -1 - outlet flow 

node number of I/O J in volume I 

node number of hot gas flow I , T 

number of H2 flows entering hot gas flow I 

po^?fon^n\TB°a?rafSon?ain?ng 9 the combustion 

p/rto hot gas now I 
node number of 02 feed J to hot gas flow 

(uncertainty estimates data input file) 

(file name variable identified m FILEIO.daij 
unused in this version 
unused in this version 
unused in this version 
• unused in this version 
■ unused in this version 

- unused in this version . * finite 

' differenc^approximatior/of "partial derivatives 

- temperature^ ractional^ncrement used in^f inite 

difference approximation of partial 
with respect to temperature 

- pressure uncertainty at node I 

- temperature uncertainty at node I 

- mass flow uncertainty at node I 

- mass flow uncertainty associated with mass fl 

- energy t uncertainty associated with voluine 1 

- uncertainty associated with volume I 


output Variable Definitions: 
No output variables 


Subroutine Calls: 

RECON 


Calling Routines: 


None 



Routine: RECON 


(reconciliation model construction and 
routing) 


Routine Function: 

Constructs a sequential quadratic programming (SQP) problem 
whose solution is the optimum reconciliation of uncertain 
test data and limited theoretical predictions for pressure, 
temperature and mass flow at specified node locations within 
the SSME flow network. Routes SQP solution logic. Outpu s 
solution of SQP problem. 


Common Blocks: 


Blank - 

TDAT - 
UDAT - 
VDAT - 
H2PRP - 

02PRP - 

H20PRP- 

TABLE - 

STD 


quadratic programming algorithm matrices 
and parameters 
test data 

uncertainty estimates 

control volume definition information 

hydrogen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 

oxygen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 

water table pressures and temperatures, 

and associated enthalpies,- entropies, and densities 

number of distinct pressures and temperatures m 

each hydrogen, oxygen, and water table .... 

reference state enthalpies entropies, and absolute 

entropies for hydrogen, oxygen, and water 


Input Variable Definitions: 

Common block inputs 


output Variable Definitions: 


I 

TTB ( IP ( I) ) 

TTB ( IT ( I ) ) 

TTB ( IW ( I ) ) 

PREC 

TREC 

WREC 

PADJ 

TADJ 

WADJ 

PPCT 

TPCT 


node number 

original test pressure at node I 

original test temperature at node I 

original test mass flow rate at node 1 

reconciled pressure at node X 

reconciled temperature at node I 

reconciled mass flow rate at node I 

pressure adjustment made at node X 

temperature adjustment made at node I 

mass flow rate adjustment made at node I 

percentage pressure adjustment made at node I 

percentage temperature adjustment made at node 

percentage mass flow rate adjustment made at 
node I 


WPCT 



Routine RECON (page 2) 

subroutine Calls: 

PROP 

CPIVOT 

Calling Routines: 


MAIN 



Routine: CPIVOT 


(solver routing routine) 


Routine Function: 

The main routing routine for the complementary pivot method, 
quadratic programming problem solver. 


Common Blocks: 

Blank - quadratic programming algorithm matrices 

and parameters 


Input Variable Definitions: 

Common block inputs 

N - dimension of square (NxN) main solver array 

N = 6*n umb er of nodes + 3*number of volumes 
+ 2*number of mass flow circuits 


Output Variable Definitions: 
Common block outputs 


Subroutine Calls: 

MATRIX 

INITIA 

NEWBAS 

SORT 


Calling Routines: 


RECON 



Routine : MATRIX 


Routine Function: 

Initializes solver inputs 


Common Blocks 
Blank 


- quadratic programming algorithm matrices 
and parameters 


Input Variable Definitions: 

Common bl °=* |£|““ on of square (NxN) main solver array 
N n ! "'number of nodes + 3*number o£ volumes 

+ ?.*number of mass flow circuits 


Output Variable Definitions: 
Common block outputs 

subroutine calls: 

None 

Calling Routines: 


CPIVOT 



Routine: INITIA 


Routine Function: 

Determines the initial almost complementary solution in the 
complementary pivot method solution strategy 


Common Blocks: 

Blank - quadratic programming algorithm matrices 

and parameters 


Input Variable Definitions: 

Common block inputs 

N - dimension of square (NxN) main solver array 

N = 6*number of nodes + 3*number of volumes 
+ 2*number of mass flow circuits 


Output Variable Definitions: 
Common block outputs 

Subroutine Calls: 

None 


calling Routines: 


CPIVOT 



Routine : NEWBAS 


Routine Function: 

Finds the new basis column to enter in terms of the current 
basis in the complementary pivot method solver 


Common Blocks: 

Blank - quadratic programming algorithm matrices 

and parameters 


Input Variable Definitions: 

Common block inputs 

N - dimension of square (NxN) main solver array 

N = 6*number of nodes + 3*number of volumes 
+ 2*number of mass flow circuits 


Output Variable Definitions: 
Common block outputs 


Subroutine Calls: 

SOLVE 

Calling Routines: 


CPIVOT 



Routine: SORT 


Routine Function: 

Finds the pivot row for the next iteration by the use of 
(simplex-type) minimum ratio rule as part of the 
complementary pivot method solver 


Common Blocks: 

Blank - quadratic programming algorithm matrices 

and parameters 

Input Variable Definitions: 

Common block inputs 

N - dimension of square (NxN) main solver array 

N = 6*number of nodes + 3*number of volumes 
+ 2*number of mass flow circuits 


Output Variable Definitions: 
Common block outputs 

Subroutine Calls: 

None 

Calling Routines: 


CPIVOT 


Routine: PIVOT 


Routine Function: 

basis^n^the^Q^vector^s^part^o^the^ompleraentar^pivot 116 

method solver 


Common Blochs: 
Blank 


quadratic programming algorithm matrices 
and parameters 


input Variable Definitions: 


Common 

N 


“^irsau :ns.;r s.s.r;r;= 

+ 2*number of mass flow circuits 


output Variable Definitions: 
Common block outputs 

Subroutine Calls : 

None 

Calling Routines: 


CPIVOT 



Routine: SOLVE 


Routine Function: 

Correlates quadratic programming problem solution as part of 
the complementary pivot method solver 

Common Blocks: 

Blank - quadratic programming algorithm matrices 

and parameters 


Input Variable Definitions: 

Common block inputs . 

m - dimension of square (NxN) main solver array 

N = 6 *number of nodes + 3*number of volumes 
+ 7 * number of mass flow circuits 


output Variable Definitions: 

Common block outputs 

Subroutine Calls: 

None 

Calling Routines: 

NEWBAS 



Routine: PROP 


Routine Function: 

Calculates hydrogen, oxygen, 


and hot gas properties 


Common Blocks: 


. ^o^table en^ies^and pities 

02PRP . denslties 

TABl£ - "umber of dis^ and water table 


input Variable 


Definitions: 


Common 

MAT 

PRSI 

TMPI 

OF 

CEFF 


block inputs 

- material type 

1 = H2, 2 = 02, 3 


hot gas 


- pressure 

- temperature 

- 02/H2 mass ratio 

- combustion efficiency 


output Variable Definitions: 


Common 

ZENTH 

ZENTR 

ZDENS 


block outputs 

- table enthalpy 

- table entropy 

- density 


Subroutine Calls: 

PRPSAT 

PRPMIX 

ITERP2 


Calling Routines: 
RECON 



Routine: PRPMIX 


Routine Function: 

Calculates hot gas mixture thermodynamic properties using a 
Dalton model 


Common Blocks: 

H2PRP - hydrogen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 
02PRP - oxygen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 
H20PRP- water table pressures and temperatures, 

and associated enthalpies, entropies, and densities 
TABLE - n umb er of distinct pressures and temperatures in 
each hydrogen, oxygen, and water table 
STD - reference state enthalpies entropies, and absolute 
entropies for hydrogen, oxygen, and water 


Input Variable Definitions: 

Common block inputs 
P - mixture pressure 

TMPI - mixture temperature 

OF - 02/H2 mass ratio 

CEFF — combustion efficiency 


Output Variable Definitions: 


Common block outputs 

HMIX - mixture enthalpy 

SMIX - mixture entropy 


PH2 

- error 

" out 

of 

range" 

PH20 

- error 

" out 

of 

range" 

P02 

- error 

" out 

of 

range" 

TMPI 

- error 

•'out 

of 

range" 


hydrogen pressure 
water pressure 
oxygen pressure 
temperature 


subroutine Calls: 


ITERP2 


Calling Routines: 


PROP 



Routine : PRPSAT 


Routine Function: 

Calculates thermodynamic properties near saturation curve 


Common Blocks: 


H2PRP - 
02PRP - 
H20PRP- 
TABLE - 
STD 


hydrogen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 

oxygen table pressures and temperatures, 

and associated enthalpies, entropies, and densities 

water table pressures and temperatures, 

and associated enthalpies, entropies, and densities 

number of distinct pressures and temperatures m 

each hydrogen, oxygen, and water table ... 

reference state enthalpies entropies, and absolute 

entropies for hydrogen, oxygen, and water 


Input Variable Definitions: 


Common block inputs 


X - pressure 

Y - temperature 

TCRT - table critical temperature . 

NX1 - n umb er of pressure values in table 

jjYl — number of temperature values in table 

NX2 

YL - table low temperature 

YH - table high temperature 

PRSl - pressure table values 

TMP1 - temperature table values 

PROP - thermodynamic property table values 

PRS2 - saturation pressure table 

TMP2 - saturation temperature table 

PROPL - saturated liquid property value 

PROPV - saturated vapor property value 


Output Variable Definitions: 

Common block outputs 

FPROP — calculated thermodynamic property 


Subroutine Calls: 


ITERP2 

ITERP1 



calling Routines 


PROP 



APPENDIX C2 

RECONCILER SOURCE CODE 
VERSION RV2-0610 



o o 


FROGRAM RECONV2 


CHARACTER * 2 4 DESC 

CHARACTER* i 2 PDAT , TDAT , UDAT , VDAT , ODAT , ODATV 


COMMON CPM( 200, 200) , CPQ( 200), LI CP, CPB( 200 , 200 ) , NL1CP, NL2CP , 

1 CPA ( 200) , NE1CP , NE2CP , IRCP, MBASIS( 300), 

2 CPW ( 200 ) , CPZ ( 200) 

C 

COMMON /VDAT/ 

1 
2 

3 

4 

5 
C 

COMMON /TDAT/ 

C 

COMMON /UDAT/ 

1 
2 
C 

COMMON /H2PRP/ 

* H2P1 (15) , H2T1 ( 11 ) 

* H2P2 ( 20 ) , H2T2 ( 11 ) 

* H2P3 ( 29 ) , H2T3 ( 25 ) 

* H2P4 ( 23 ) , H2T4 ( 25 ) 

COMMON /02PRP/ 

* 02P1(13) , 02T1 ( 1 6 ) 

* 02P2(13) , 02T2 ( 17 ) 

* 02P3 ( 5 ) , 02T3 ( 61 ) 

COMMON /H20PRP/ 

* H20P1 ( 7 ) , H20T1 ( 13 
C 

COMMON /TABLE/ 

* NH2P( 4 ) ,NH2T(4) ,N02P(3) ,N02T(3) , NH20P( 1 ) ,NH2OT( 1 ) 

COMMON /STD/ 

* HH2REF , H02REF , HWAREF , SH2REF , S02REF , SWAREF , SH2A , S02A , 

* SWAA 
C 

DIMENSION 

* NH2PA( 4 ) ,NH2TA( 4 ) ,N02PA( 3 ) ,N02TA( 3 ) ,NH20PA( 1 ) , NH20TA( 1 ) 

C 

CHARACTER* 70 PTITLE 
C 

DATA ( NH2PA( I ) , 1-1 , 4 ) /15 , 20 , 29 , 23/ 

DATA ( NH2TA( J) , J-l , 4 )/ll , 11 , 25 , 25/ 

DATA (N02PA(I),I-1,3)/13,13,5/ 

DATA (N02TA(J) , J-l , 3 )/16 , 17 , 61/ 

DATA (NH20PA(I) ,1-1,1) /7/ 

DATA ( NH20TA( J) ,J— 1,1)/13/ 

C 

DO 90 1-1,4 
NH2P( I ) -NH2PA( I ) 

90 NH2T( I ) -NH2TA( I ) 

DO 91 1-1,3 
N02P( I ) — N02PA( I ) 

91 N02T( I )— N02TA( I ) 

NH20P( 1 )=NH20PA( 1 ) 


, H2H1 (15,11) ,H2S1(15,11) , H2D1 (15,11) , 
, H2H2 (20,11) , H2S2 (20,11) ,H2D2( 20,11) , 
, H2H3 (29,25) , H2S3 ( 29 , 25 ) ,H2D3 ( 29 , 25 ) , 
, H2H4 (23,25), H2S4 (23,25) ,H2D4 (23,25) 

, 02H1( 13, 16) ,0251(13,16) ,0201(13,16) , 
, 02H2( 13,17) ,02S2( 13, 17) ,0202(13,17) , 
, 02H3 (5,61), 02S3 (5,61), 02D3(5,61) 

) ,H20H1(7,13) ,H20S1(7,13) ,H20D1(7,13) 


IENV, IPCTTH , MAXSTG , NHG , NTMFC, NTNOD, NTVOL , 
IA( 20 ) , IP( 20 ) , IT( 20 ) , IW( 20 ) , MAT(20), 

MIO ( 5 ) , MODIR( 5,20) , IMFCN(5,20), 

NIO( 5 ) , IODIR (5,20) , IVOLN(5,20), 

NH2HG( 5 ) , N02HG ( 5 ) , NODHG ( 5) , ICEFF(5), 

IH2HG( 5,5), I02HG( 5,5) 

TTB(IOO), NDESC , NTTB, DESC(5) 

IPRPD , ITTBD , IROWA , JCOLA , ITORDX, ITORDY , 
DPF , DTF , UP( 20 ) , UT(20), UW(20), 

UEVOL( 5 ) , USVOL( 5 ) , UWMFC(5) 



o o o 


NH20T ( 1 )=NH20TA( 1 ) 


HH2REF - 
H02REF = 
HWAREF * 
SH2REF = 
S02REF = 
SWAREF = 
SH2A 
S02A 
SWAA 


1790 . 091 
234 . 681 
1339.990 
15.440 
1.530 
2 
15 
1 
0 


294 

481 

531 

928 


OPEN 

READ 

OPEN 

OPEN 

OPEN 

OPEN 


( 


FILE 
* ) 


- 'FILEIO.DAT' , STATUS - 'OLD' ) 
PDAT, TDAT , UDAT, VDAT , ODAT, ODATV 


C 

c 

c 


8, FILE 

12, FILE 

13, FILE 

14, FILE 


PDAT 

TDAT 

UDAT 

VDAT 


STATUS 

STATUS 

STATUS 

STATUS 


'OLD' 

'OLD' 

'OLD' 

'OLD' 


OPEN ( 21, FILE 
OPEN ( 22, FILE 


ODAT ) 
ODATV ) 


** 


READ IN H2 PROPERTY TABLE INTO ARRAYS 


** 


DO 10 ITBL*1 , 4 

READ( 8 , 902 ) PTITLE 
DO 10 1=1 ,NH2P( ITBL) 

DO 10 J=1 ,NH2T( ITBL) 

IF( ITBL.EQ. 1) READ ( 8 , * ) H2P1(I) ,H2T1( J) , 

1 H2H1(I,J) ,H2S1(I,J) ji 

IF( ITBL.EQ. 2) READ ( 8 , * ) H2P2( I ) ,H2T2 ( J) , 

, H2H2 f I J) H2S2 ( I # J) t H2D2 ( I , J) 

IF (ITBL.EQ. 3) READ( 8 , * ) H2P3 ( I ) ,H2T3 ( J) , 

1 H2H3 1 J) H2S3(I,J),H2D3(I,J) 

1 IF ( ITBL. EQ. 4 ) READ ( 8 , * ) «2P4( I) ,H2T4( J) , 

1 H2H4(I,J) ,H2S4(I,J) ,H2D4(I,J) 

10 CONTINUE 

** READ IN 02 PROPERTY TABLE INTO ARRAYS ** 

DO 20 ITBL*1,3 

READ(8,902) PTITLE 
DO 20 1=1 ,N02P( ITBL) 

DO 20 J=1,N02T( ITBL) 

IF(ITBL.EQ.l) READ ( 8 , * ) 02P1(I) ,02T1( J) . 

1 02H1(I,J) ,02S1(I,J) ,02D1(I»J) 

IF (ITBL.EQ. 2) READ ( 8 , * ) 02P2(I) ,02T2( J) , 
i n2H2fI J ) 02S2 ( I , J 1 t 02D2 ( I f J) 

IF( ITSL. EQ . 3 ) READ ( 8 , * ) 02P3 ( I ) ,02T3 ( J) #■ 

1 02H3 (I,J) ,02S3(I,J) ,02D3 ( I , J) 

20 CONTINUE 

** READ IN STEAM PROPERTY TABLES INTO ARRAYS ** 


C 

C 



c 


DO 30 ITBL =1,1 


READ (8, 90 21 PTITLE 

DO 30 1=1, NH20P ( ITBL ) 

DO 30 J = 1, NH20T( ITBL) 

IF ( ITBL . EQ. 1 ) READ ( B , * ) H20P1 ( I ) , H20T1 ( J) 
1 H20H1(I,J) ,H20S1(I,J) , H20D1 ( I , J ) 


30 CONTINUE 


READ 

READ 

READ 


( 12 ,*) 

( 12 ,*) 

( 12 ,*) 


NDESC, NTTB 
( DESC( I ), 
( TTB( I ) , 


WRITE (21,901) ( DESC( I ), 1 


NDESC 

NTTB 

NDESC 


READ 

READ 

2 

3 

4 

5 

6 
7 


(14,*) 

(14,*) 


IENV , 


IPCTTH , MAXSTG, NHG , NTMFC, NTNOD , NTVOL 


( 


IA( I ) , 
IP(I) , 
IT( I) , 
IW(I) , 
MAT(I) , 
MIO(I) , 
NIO(I) , 


I 

I 

I 

I 

I 

I 

I 


DO 50 1 = 1 

READ (14,*) 

1 

50 CONTINUE 


NTMFC 

( MODIR( I , J) , 
( IHFCN( I , J) , 


60 


DO 60 1=1, NTVOL 

READ (14,*) ( IODIR(I,J), 

! ( IVOLN ( I , J) , 

CONTINUE 


1 

1, 

1, 

1, 

1, 

1, 

1 , 


J 

J 


J 

J 


NTNOD 

NTNOD 

NTNOD 

NTNOD 

NTNOD 

NTMFC 

NTVOL 


), 
), 
), 
) , 
) , 
) , 
) 


1, 

1, 


1, 

1, 


MIO(I) 

MIO(I) 


NIO(I) 

NIO(I) 


), 

) 


l. 

) 


IF ( 
READ 


NHG . GT. 0 ) THEN 
(14,*) ( NODHG(I) , I “ 

( NH2HG ( I ) , I ■ 
( N02HG( I ) , I = 
( ICEFF(I), I = 
- 1, NHG 

*) ( IH2HG( I , J) , 

( I02HG( I , J) , 


70 


DO 70 I 
READ (14, 
L 

CONTINUE 

ENDIF 


1, 

1, 

1, 

1, 

J 

J 


NHG 

NHG 

NHG 

NHG 

= 1 , 

= 1. 


), 

), 

), 

) 

NH2HG ( I ) 
N02HG ( I ) 


), 

) 


READ (13,*) IPRPD, 
DPF, DTF, 


ITTBD, IROWA, JCOLA, ITORDX, ITORDY, 


UP(I) , 

UT( I ) , 
UW(I) , 
UWMFC(I) , 
UEVOL(I) , 
USVOL(I) 


I 

I 

I 

I 

I 

I 


1, 

1. 

1, 

1, 

1, 

1, 


NTNOD 

NTNOD 

NTNOD 

NTMFC 

NTVOL 

NTVOL 


CALL RECON 
*901 FORMAT ( 


10 ( /, IX, A24 ) ) 



902 FORMAT ( /A70/ ) 
C 
C 

STOP 

END 

C* *.•*»*«.»*«********* 

SUBROUTINE RECON 


C 

c 

c 


c 


c 


c 


c 

c 


c 


c 


c 

c 

c 

c 


RECON - RECONCILIATION 


CHARACTER* 2 4 DESC 
REAL JOULE 


DIMENSION 

1 

2 

3 

4 


DDDPN ( 20 ) , DDDTN ( 20 ) 
DSDPN( 20 ) , DSDTN ( 20 ) 
ASTD( 20 ) , HSTD( 20 ) 
REVA( 20 ) , REVP( 20 ) 

CPQQT( 100 , 100 ) 


, DHDPN( 20) , DHDTN( 20 ) 
, DENS (20) , 

, SSTD( 20) , 

, REVT( 20 ) , REVW( 20 ) 
TTBREV( 100) 


A(100 ,100) 


COMMON CPM( 200 , 200 ) , CPQ( 200), 
CPA( 200 ) , NE1CP , NE2CP , 
i CPW( 200 ) , CPZ ( 200 ) 


L1CP , CPB( 200 , 200 ) , NL1CP , 
IRCP , MBASIS( 300) , 


NL2CP , 


COMMON /VDAT / 

1 

2 

3 

4 

5 


IEHV, IPCTTH, HAXSTG, NHG , NTMFC, NTNOD, 
IA( 20 ) , IP( 20 ) , IT( 20 ) , IW( 20 ) , MAT(20), 
MIO( 5) , MODIR( 5,20) , IHFCN(5,20), 

NIO( 5 ) , IODIR( 5,20), IVOLN(5,20) 

NH2HG( 5) , N02HG( 5 ) , NODHG( 5) , 

IH2HG( 5,5), I02HG( 5,5) 


ICEFF(5) 


NTVOL, 


COMMON /TDAT/ TTB(IOO), NDESC , NTTB , DESC(5) 


COMMON /UDAT/ IPRPD, ITTBD , IROWA , JCOLA , 
DPF , DTF , UP( 20 ) , UT( 20) , 
; UEVOL( 5 ) , USVOL( 5 ) , UWMFC(5) 


ITORDX, ITORDY, 
UW( 20) , 


COMMON /H2PRP/ 

1 H2P1(15) ,H2T1(11) 

2 H2P2 ( 20 ) , H2T2 ( 11 ) 

3 H2P3 ( 29 ) , H2T3 ( 25 ) 

4 H2P4 ( 23 ) ,H2T4 ( 25 ) 
COMMON /02PRP/ 

1 02P1 ( 13 ) , 02T1 ( 16 ) 

2 02P2 ( 13 ) , 02T2 ( 17 ) 

3 02P3 ( 5 ) , 02T3 ( 61 ) 
COMMON /H20PRP/ 

1 H20P1(7),H20T1(13 


, H2H1 ( 15 
,H2H2( 20 
,H2H3 ( 29 
, H2H4 ( 2 3 

,02H1( 13 
,02H2( 13 
, 02H3 ( 5 , 


,11) ,H2S1 ( 15 
,11) , H2S2 (20 
,25) ,H2S3( 29 
,25) ,H2S4 ( 23 

,16) ,02S1( 13 
,17) ,02S2(13 
61), 02S3 ( 5 , 


,11) ,H2D1(15,11) , 
,11) ,H2D2(20,11) , 
,25) ,H2D3(29,25) , 
,25) ,H2D4 (23,25) 

,16) ,02D1(13,16) , 
,17) ,0202( 13 ,17) , 
61), 02D3(5,61) 


) , H20H1 (7,13) ,H20S1 (7,13) ,H20D1(7,13) 


1 C HH2REF^H02REF , HWAREF , SH2REF , S02REF , SWAREF , SH2A, S02A , 
2 SWAA 

1^NH2P( 4^ , NH2T^ 4 ) ,N02P(3) ,N02T(3) ,NH20P(1) ,NH20T(1) 
PARAMETER ( JOULE - 778.16, GC * 32.174 ) 


ISTG - 1 

TENV = TTB( IENV ) 


DO 10 I 
REVMI) 
REVP ( I ) 
REVT( I ) 
REVW ( I ) 
CONTINUE 


1, NTNOD 
TTB( IA ( I ) ) 
TTB ( IP(I) ) 
TTB( IT ( I ) ) 
TTB( IW(I) ) 


DO 70 1=1, NTNOD 

p = REVP ( I ) 

T = REVT( I ) 

W = REVW { I ) 


IF ( MAT(I) 
IF ( MAT(I) 

CALL PROP ( 
DENS ( I ) = 

HN 

HSTD(I) = 
SN 

SSTD(I) * 
AN 

ASTD(I) “ 


.GE. 3 ) GO TO 40 
.GE. 2 ) GO TO 30 

1, P, T, 0.0, 0.0, H, S, 
RHO 

H - HH2REF 
HN 

S - SH2REF + SH2A 
SN 

HN - TENV * SN 
AN 


RHO) 


P2 

CALL PROP ( 
HN2 

SH2 * 

DDDPN(I) = 
DHDPN(I) - 
DSDPN(I) - 


p + DPF * P 

1 t p2 , T , 0.0, 0.0, H2, S2, RH02) 
H2 - HH2REF 
S2 - SH2REF + SH2A 
( RHO 2 - RHO ) / ( P2 - P ) 

( HN2 - HN ) / ( P2 - P ) 

( SN2 - SN ) / ( P2 - P ) 


T2 

CALL PROP 

HN2 

SN2 

DDDTN(I) 

DHDTN(I) 

DSDTN(I) 


T + DTP * T 

1, P, T2 , 0.0, 0.0, H2 , S2 , 
H2 - HH2REF 
S2 - SH2REF + SH2A 
( RHO 2 - RHO ) / ( T2 - T ) 

( HN2 - HN ) / ( T2 - T ) 

( SN2 - SN ) / ( T2 - T ) 


RHO 2 ) 


30 


GO TO 70 


CALL PROP ( 
DENS(I) - 
HN 

HSTD(I) - 
SN 

SSTD(I) ” 

AN 

ASTD(I) = 


2, P, T, 0.0, 0.0, H, S, RHO) 
RHO 

H - H02REF 
HN 

S - S02REF + SO 2 A 
SN 

HN - TENV * SN 
AN 


P2 

CALL PROP ( 
HN2 

SN2 — 

DDDPN ( I ) - 

DHDPN(I) - 
DSDPN(I) = 


P + DPF * P 

2, P2 , T, 0.0, 0.0, H2 , S2. RH02) 
H2 - H02REF 
S2 - S02REF + SO 2 A 
( RHO 2 - RHO ) / ( P2 - P ) 

( HN2 - HN ) / ( P2 - P ) 

( SN2 - SN ) / ( P2 - P ) 



CALL PROP ( 2. P. T2 ,0.0, 0.0, H2 , S2, RH02) 

£’ I g I S02REF ♦ S02A _ 

DDDTN(I) * ( RH02 - RHO ) / _ 

DHDTN(I) = ( HN2 - HN / < T2 i i 

DSDTN(I) = ( SN2 _ SN ' 


/ ( T2 - T ) 


GO TO 70 

40 DO 42 IDHG = 1, NHG THEN 

IF ( NODHG(IDHG) .EQ. I ) THEN 
IH G - IDHG 
GO TO 44 
ELSE 
ENDIF 

42 CONTINUE 

To\o = IH2iS - 1. KH2HG( IHG ) 

g“ : S??l«55< JnUM 

CONTINUE 

ur>5 at 0 » 0 

DO 60 I02IN - 1. } 

NNUM = I02HG( IHG, ^ 

„02 = W02 + REVW( NNUM ) 

CONTINUE 


44 


50 


60 


TTB( ICEFF( IHG ) ) 

<°V, SfoF. CEFF, HMIX, SMIX, DMIX) 

DMIX* 

HMIX 

SMIX 

HMIX “ TENV * SMIX 
AMIX 


CEFF 

OF “ 

CALL PROP ( 

DENS ( I ) - 

HSTD(I) - 
SSTD(I) - 
AMIX 

ASTD(I) - 

P2 oonp T 4 P U DP t *OF, CEFF, H2MIX, S2MIX, D2MIX) 

CALL PROP ( 4 , P2 , 1 , ur * ^ . 

DDDPN(I) • ( D2MIX - DHIX W f P2 - P ) 

DHDPN ( 1 1 - ( H2MIX - HMIX / ” . p ) 

DSDPN(I) ■ ( S2HIX SHXA ) / \ 

T2 oort p 7 4 T P °T2 # OF, CEFF, H2MIX, S2MIX, D2MIX) 

CALL PROP ( 4, P, T2, ur, > __ _ T \ 

DDDTN(I) - ( D2MIX - DMIX ! ( t2 - T ) 

DHDTN(I) - ( H2ME * ” Ztx W ( T2 - T ) 

DSDTN(I) - ( S2MIX - SMIX ) / ( T2 


70 CONTINUE 


M 

N 


3 * NTNOD + 2 * NTMFC + 3 * NTVOL 
3 * NTNOD 


80 


DO 80 I • 1# M 
DO 80 J “ N 
A(I,J) “ °*° 

CONTINUE 


C 



DO 8 2 I = 1 , N 
DO 8 2 J = 1 , N 
CPQQT ( I , J ) ■ 0.0 

82 CONTINUE 


DO 84 ITNOD 

11 

12 
13 

CPQQT ( II .ID 
CPQQT (12,12) 
CPQQT ( 13,13) 
CPQ( II ) 
CPQ(I2) 

CPQ( 13 ) 

84 CONTINUE 


, NTNOD 
ITNOD 

ITNOD + NTNOD 
ITNOD + 2 * NTNOD 
4 . / UW( ITNOD ) ** 2 
4. / UP( ITNOD ) ** 2 
4 . / UT( ITNOD ) ** 2 

-4. / UW( ITNOD ) 

-4. / UP( ITNOD ) 

-4. / UT( ITNOD ) 


DO 90 ITNOD 

11 

12 
13 

Mil, ID 

M 12 , 12 ) 

A( 13 , 13 ) 
CPQ(N+I1 ) * 

CPQ( N+I2 ) a 
CPQ( N+I3 ) * 

90 CONTINUE 


- 1, NTNOD 
ITNOD 

ITNOD + NTNOD 
ITNOD + 2 * NTNOD 
- 1 . 

- 1 . 

- 1 . 

2. * UW( XI ) 

2. * UP( II ) 

2. * UT( II ) 


C 

c 


NQROW 

= 2 * N 

DO 98 

ITMFC ® 1 0 NTMFC 

SUMM 

- 0. 

XI 

- ITMFC 

12 

« ITMFC + NTMFC 


DO 92 
IOD 
IMN 
W 

SUMM 


110 - 1, MIO( ITMFC ) 

» HODIR( ITMPC , IIO ) 

- IHFCN ( ITMFC, IIO ) 

- REVW( IMN ) 

» SUMM + IOD * ( W - UW( IMN 


92 CONTINUE 


) 


C 


c 


94 


C 


98 


C 


DO 94 IIO * 1» 

IOD 

IMN 

A(N+I1,IKN) - 
A( N+I2 , IMN) - 
CONTINUE 

CPQ( NQROW+I1 ) 
CPQ( NQROW+I2 ) 
CONTINUE 


IO( ITMFC ) 

MODIR ( ITMFC, IIO ) 
IHFCN ( ITMFC, IIO ) 
IOD 


. SUMM + UWMFC ( ITMFC ) 
a - SUMM + UWMFC ( ITMFC ) 


NAROW 

NQROW 


N + 2 * NTMFC 
2 * N + 2 * NTMFC 


DO 110 

11 

12 


ITVOL “ 1 , NTVOL 
. ITVOL 

m ITVOL + NTVOL 



I3 = ITVOL + 2 * KTVOL 

CEO * 0 • 

CSO = 0 • 

SUMMI = 0 • 


DO 99 

IOD 

IVN 

AREA 

W 

RHO 

EHTH 

ENTR 

AVAIL 


XIO - 1 , HIO( ITVOL ) 
= IODIR( ITVOL, 1 10 
= IVOLN( ITVOL, IIO 
= REVA( IVN ) 

= REVW ( IVN ) 

= DENS ( IVN ) 

= HSTD ( IVN ) 

= SSTD( IVN ) 

= ASTD( IVN ) 


) 

) 


99 


£0i o IOD * W * ENTH 
f\ AREA .GT. 0. ) THEN 

* v t 1 / 


ELSE 

CE02 

ENDIF 

CEO 

CSO 

CONTINUE 


- 0 . 


CEO + 
CSO + 


CE01 + CE02 
XOD * W * ENTR 


JOULE * RHO** 2 


C 


C 


SUMEQ - 0 • 0 

SUMSQ = 0 * 0 


DO 100 

IOD 

IVN 

AREA 

W 

P 

T 

RHO 

ENTH 

ENTR 

AVAIL 

DDDP 

DHDP 

DSDP 

DDDT 

DHDT 

DSDT 


IIO - lr 
3 IODIR( 
. IVOLN ( 


NI0( ITVOL 
ITVOL, IIO 
ITVOL, IIO 


- REVA( IVN ) 
. REVW( IVN ) 

- REVP( IVN ) 
. REVT( IVN ) 

- DENS( IVN ) 
=, HSTD( IVN ) 
O SSTD( IVN ) 


» ASTD( IVN ) 

. 0DDPN( IVN ) 
3 DHDPN( IVN ) 
= DSDPN( IVN ) 
=. DDDTN( IVN ) 
» DHDTN ( IVN ) 
= DSDTN( IVN ) 


) 

) 

) 


CEM 

CEP 

CET 

CSM 

CSP 

CST 


IOD * EHTH 
IOD * W * DHDP 
IOD * W * DHDT 
IOD * EHTR 
IOD * W * DSDP 
IOD * W * DSDT 


IF 


1 

1 

1 


( AREA 
CEM * 

CEP - 

CET - 


T • 0 . ) THEN 

“ VkV “ o ^“"» k >** 2 * “ ,EA " ! > 

rp - IOD * W**3 * DDDP / _ N 

GC * JOULE * RHO**3 * ***** 2 > 

:ET - IOD * W**3 * DDDT / 

GC * JOULE * RHO** 3 * AREA** 2 ) 


AREA** 2 ) 



ENDIF 


C 


C 


c 


NR0W1 

a(nroki,:vn) 

A( NROW1 , IVN+NTNOD) 
MNROW1 ,IVN+2*NTNOD) 


NAROW +11 

CEM 

CEP 

CET 


NROW2 

A( NROK2 , IVN ) 

A( NROW2 , IVN+NTNOD) 

A( HROK2 , IVN+ 2 * NTNOD ) 


= NAROW +12 
= -CEM 

- -CEP 

- -CET 


NROW3 

A( NROW3 , IVN) 

A( NROW3 , IVN+NTNOD) 
A(NROW3 , IVN+2+NTNOD) 


NAROW+I3 

-CSM 

-CSP 

-CST 


SUMEQ = SUMEQ + CEM * UW( IVN ) + CEP * UP( IVN ) + 
CET * UT( IVN ) 

SUMSQ = SUMSQ + CSM * UW( IVN ) + CSP * UP( IVN ) + 
, CST * UT( IVN ) 


C 

100 CONTINUE 


C 


C 

110 

C 


120 

C 


130 

c 


140 


CPQ ( NQROW+I 1 ) 
CPQ ( NQROW+I 2 ) 
CPQ( NQROW+I 3 ) 


CEO + UEVOL( ITVOL ) - SUMEQ 

- CEO + UEVOL( ITVOL ) + SUMEQ 

- CSO + USVOL( ITVOL ) + SUMSQ 


CONTINUE 

DO 120 I - 1, N 
DO 120 J » 1, N 
CPM(I,J) “ CPQQT ( I, J ) 
CONTINUE 

DO 130 I * 1 ( H 
DO 130 J - 1, N 
CPM(N+I,J) * A( I, J ) 
CPM(J,N+I) » - A( I, J ) 

CONTINUE 


NPl a N + 1 

MN = M + N 

DO 140 I - NPl, MN 

DO 140 J » NPl, MN 

CPM( I , J) “ 0.0 

CONTINUE 


C 


C 


C 


300 


CALL CPIVOT ( HN ) 
IF ( IRCP .EQ. -1 ) 
IF ( IRCP .EQ. -2 ) 
IF ( IRCP .EQ. -3 ) 


WRITE ( 21, 986 ) 
WRITE ( 21, 987 ) 
WRITE ( 21, 988 ) 


WRITE ( 21, 981 ) 
WRITE ( 21, 982 ) 


DO 300 I ■ 1, NTNOD . 

WRITE ( 21, 951 ) I, TTB(IP(I)), TTB(IT(I)), 

CONTINUE 


TTB ( IW ( I ) ) 


C 



310 


WRITE 

( 

21, 983 ) 



WRITE 

( 

21, 982 ) 



DO 210 


1=1, NTNOD 

) + ( 

REVW ( 

I 

) = REVW ( 

I 

REVP ( 

T 

) = REVP ( 

I 

) + ( 

REVT( 

T 

) = REVT ( 

I 

) + ( 

WREC 

= 

REVW ( I ) 



PREC 

= 

REVP ( I ) 



TREC 

= 

REVT( I ) 


PREC, 

WRITE 

( 

21, 951 ) 

I, 

CONTINUE 



WRITE 

( 

21, 984 ) 



WRITE 

( 

21, 982 ) 



DO 320 

1=1, NTNOD 

TTB( 3 

WADJ 

s 

REVW( I ) 

— 

PADJ 

= 

REVP( I ) 

“ 

TTB ( 3 

TADJ 


REVT( I ) 

- 

TTB( 3 

WRITE 

( 

21, 951 ) 

I ( 

, PADJ 


) - UW( I ) ) 

+ ntnod ) - UP( I ) ) 

+ 2 * NTNOD ) “ UT( I ) ) 


I ) ) 


320 CONTINUE 


WRITE ( 21, 985 ) 
WRITE ( 21, 982 ) 


C 


330 

C 


DO 3 3 0 1 * 1, NTNOD 

WPCT = 100 * ( REVW( I ) 

PPCT * 100 ^ ( REVP( I ) 

TPCT - 100 * ( REVT ( I ) 

WRITE ( 21, 951 ) I, PPCT 
CONTINUE 


- TTB( IW( I ) 

- TTB( IP( I ) 

- TTB( IT( I ) 
, TPCT, WPCT 


) ) / TTB( IW( 
) ) / TTB( IP( 
) ) / TTB( IT( 


I ) 
I ) 
I ) 


) 

) 

) 


IF ( ISTG .LT. MAXSTG ) THEN 
ISTG * ISTG + 1 
GO TO 20 
ELSE 
ENDIF 


500 CONTINUE 

TTBREV( IENV ) 
TTBREV ( IPCTTH ) 


TTB( IENV ) 
TTB( IPCTTH ) 


C 


510 

C 


520 

C 


DO 510 1*1, NTNOD 

TTBREV ( IA( I ) ) * 

TTBREV ( IP(I) ) * 

TTBREV ( IT(I) ) * 

TTBREV ( IW(I) ) “ 


REVA( I ) 
REVP( I ) 
REVTt I ) 
REVW( I ) 


CONTINUE 


DO 520 1*1, NHG , . 
TTBREV ( ICEFF(I) ) ■ TTB( ICEFF(I) ) 
CONTINUE 


WRITE ( 22, * ) 
WRITE ( 22, 990 ) 
WRITE ( 22, 989 ) 


NDESC, NTTB 
( DESC(I) , 
( TTBREV ( I ) , 


I - 1, NDESC ) 
I - 1, NTTB ) 


C 

951 


FORMAT ( 9X , 16, 3F15.2 ) 



981 FORMAT ( 

982 FORMAT ( 

983 FORMAT ( 

984 FORMAT ( 

985 FORMAT ( 

986 FORMAT ( 

987 FORMAT ( 

988 FORMAT ( 

989 FORMAT ( 

990 FORMAT ( 
C 

RETURN 


,/:X, 'RECONCILED NODE DATA ) 

/ / 2 X , 'BALANCING ADJUSTMENTS ) 

: P N ?Sp,.rSST™’‘sO I .iTIOH OBTAINED ' ) 

'//lx', 'TRIVIAL COMPLEMENTARY SOLUTION' ) 
//2X, 'SOLUTION OBTAINED' ) 

6F12.2 ) 

2X , 1H ' , A24 , 1H ' ) 


SUBROUTINE CPIVOT (N) 


COKHON AM(200, 200) ,0(200), LI, B(200 200) ,»L1.HL2,A( 200) ,Htl,NE2. 

1 IR,HBASIS( 300) , W( 200 ) , Z ( 200 ) 


description ^^^ N aSTy ON containing the 

q ffSS SUBSCRIPTED ARRAY CONTAINING THE 

Ll S E INTEGER F VARIABLE 0 INDICATING THE NUMBER OF 

L1 ITERATIONS TAKEN FOR EACHPROBLEM. 

B A TWO DIMENSIONAL ARRAY CONTAINING THE 

w KSB %&££%& 

z si NGLY^ SUBSCRIPTED 0 ARRAy'^CONTAINING THE VALUES 

2 OF Z VARIABLES IN EACH SOLUTION. 

mti XN INTEGER VARIABLE TAKING VALUE 1 OR 2 DEPEND 

KL1 i^STSSlHE* VARIABLE W OR 2 LEAVES THEBAMS 

SIMILAR TO NL1 BUT INDICATES VARIABLE QfTERim* 
Sl 2 jS INTEGER VARIABLE INDICATING WHAT COMPONENT 

OF W OR Z VARIABLE LEAVES THE BASIS. 

SlHllS. TO NL2 BUT INDICATES VARIABLE ENTERING 
A A SINGLY SUBSCRIPTED ARRAY CONTAINING 1 THE 

ELEMENTS OF THE TRANSFPORMED COLUMN THAT IS 

IR S T I»TE=™«S^E DEMOTING THE PIVOT ROM AT 

EACH ITERATION- ALSO USED TO INDICATE 
ALGORITHM TERMINATION 

IR - -2 COMPLEMENTARY SOLUTION DETERMINED 
IR - -1 PROBLEM HAS NO COMPLEMENTARY SOLUTION 

MHXSIS A SINGLY SUBSCRIPTED ARRAY -INDICATOR FOR THE 
MBASIS y^j^LES. TWO INDICATORS ARE USED FOR 

EACH BASIC VARIABLE-ONE INDICATING WHETHER 

^ ?s f £ orTand another indicating what 

COMPONENT OF W OR Z. 

PROGRAM CALLING SEQUENCE, N IS THE SIZE OF MATRIX AM 


CALL MATRIX ( N ) 

CALL INITIA ( N ) 

IF ( IR • EQ. -2 ) GO TO 5 


oo OOOO ODD 


c 

4 CALL NEWBAS ( N ) 

IF ( IR .EQ. -3 ) GO TO 5 
C 

CALL SORT ( N ) 

IF ( IR -EQ- "I ) GO TO 5 
C 

CALL PIVOT ( N ) 

GO TO 4 
C 

5 RETURN 
END 

C 

C 

C 

SUBROUTINE MATRIX ( N ) 


C PURPOSE - TO INITIALIZE THE VARIOUS INPUT DATA 


C 


COMMON AM( 200 , 200 ) , Q(200), LI, B( 200, 200) 
NE1, NE2, IR, MBASIS ( 300 ) , W(200), 


, NL1, 
Z( 200 ) 


NL2, 


C IN ITERATION 1, BASIS INVERSE IS AN IDENTITY MATRIX 

C 

DO 5 J » 1, N 

DO 4 I ■ 1, N 

IF ( I .EQ. J ) GO TO 3 

B(I,J) - 0.0 

GO TO 4 

3 B(I,J)-1.0 

4 CONTINUE 

5 CONTINUE 
C 

RETURN 

END 


A( 200 ) , 


SUBROUTINE INITIA ( N ) 

PURPOSE - TO FIND THE INITIAL ALMOST COMPLEMENTARY SOLUTION 
BY ADDING AN ARTIFICIAL VARIABLE ZO. 


COMMON AM( 200 , 200 ) , Q(200), LI, B(200,200) 
NE1, NE2, IR, MBASIS ( 300) , W(200), 


, NL1, 
Z ( 200 ) 


NL2, 


A( 200 ) , 


SET ZO EQUAL TO THE MOST NEGATIVE Q(I) 

I - 1 

J - 2 

1 IF ( Q( I ) .LE. Q( J) ) GO TO 2 
I - J 

2 J = J + 1 

IF ( J -LE. N ) GO TO 1 

UPDATE Q VECTOR 
IR - I 

T1 “ -Q( IR ) 

IF ( T1 .LE. 0.0 ) GO TO 9 

DO 3 I * 1, N 
Q ( I ) - Q( I ) + T1 



on o o ooo noon non 


3 CONTINUE 

Q(IR) = T1 
C 

C UPDATE BASIS INVERSE AND INDICATOR VECTOR 
C OF BASIC VARIABLES 


DO 4 J * 1 , N 

B(J,IR) 

-1. 0 

W( J) 

Q( J ) 

2(J) 

0.0 

MBASIS(J) = 

1 

L 

N + J 

MBASIS(L) = 

J 

CONTINUE 


NL1 

1 

L 

N + IR 

NL2 

IR 

MBAS IS (IR) = 

3 

MBASIS(L) 

0 

W(IR) 

0.0 

20 

Q( IR ) 

LI 

l 

RETURN 


IR * -2 


RETURN 


END 



SUBROUTINE NEWBAS ( N ) 

PURPOSE - TO FIND THE NEW BASIS COLUMN TO ENTER IN 
TERMS OF THE CURRENT BASIS. 

COMMON AM( 200 , 200 ) , Q(200), LI, B(200,200), NL1, NL2 , A(200), 
1 NE1, NE2, IR, MBASIS( 300 ) , W(200), Z(200) 

IF NL1 IS NEITHER 1 NOR 2 THEN THE VARIABLE Z0 LEAVES THE 
BASIS INDICATING TERMINATION WITH A COMPLEMENTARY SOLUTION 
IF ( NL1 .EQ. 1 ) GO TO 2 
IF ( NL1 .EQ. 2 ) GO TO 5 

CALL SOLVE ( N ) 

IR - -3 

RETURN 

2 NE1 = 2 

NE2 = NL2 

UPDATE NEW BASIC COLUMN BY MULTIPLYING BY BASIS INVERSE. 

DO 4 I • 1, N 

T1 * 0.0 

DO 3 J - 1, N 

3 T1 - T1 - B( I, J ) * AM( J, NE2 ) 

A( I ) = T1 

4 CONTINUE 
RETURN 

C 



o o o o non 


5 NE1 = 1 

NE2 = NL2 

DO 6 I * 1 . N 

A( I ) = B ( I, NE2 ) 

6 CONTINUE 

RETURN 

END 


SUBROUTINE SORT ( N ) 


PURPOSE - TO FIND THE PIVOT ROW 
USE OF (SIMPLEX-TYPE) 


FOR THE NEXT ITERATION BY 
MINIMUM RATIO RULE. 


COMMON AM( 200 , 200 ) , Q(200), 
NE1, NE2, IR, MBAS IS ( 


LI, B( 200, 200), NL1 , 
300) , W ( 200 ) , Z ( 200 ) 


NL2, 


1=1 

1 IF ( A( I ) -GT. 0.0 ) GO TO 2 
I =1 + 1 

IF ( I .GT. N ) GO TO 6 
GO TO 1 


C 

2 T1 = Q( I ) / A( I ) 

IR = I 

3 1 =1 + 1 

IF ( I -GT. N ) GO TO 5 

IF ( A( I ) -GT. 0.0 ) GO TO 4 

GO TO 3 


4 T2 = Q( I ) / M 1 ) 

IF ( T2 .GE. T1 ) GO TO 3 
IR = I 

T1 = T2 

GO TO 3 


C 

5 RETURN 


C FAILURE OF THE RATIO RULE INDICATES TERMINATION WITH 
C NO COMPLEMENTARY SOLUTION. 

C 

6 IR = -1 ' 

RETURN 

END 

C 

C 

C 

SUBROUTINE PIVOT ( N ) 


C PURPOSE - TO PERFORM THE PIVOT OPERATION BY UPDATING THE 
C THE INVERSE OF THE BASIS AND Q VECTOR. 

° COMMON AM( 200, 200) , Q(200), LI, 8(200,200), NLX, NL2 

1 NE1, NE2 , IR, MBASIS ( 300 ) , W(200), Z(200) 


DO 1 I = 1, N 

1 B( IR , I ) * B( IR, I ) / A( IR ) 
Q( IR) = Q( IR ) / A( IR ) 


THE 

A ( 200 ) , 


, A( 200 ) , 



non o o o 


DO 3 I = 1 , N 

IF ( I .EQ. IR ) GO TO 3 

Q ( I ) = Q( I ) - Q( IR ) 


A ( I ) 


DO 2 J = 
B(I,J) = 

2 CONTINUE 

3 CONTINUE 


1, N 

B( I, J ) 


B( IRi J ) * 


A ( I ) 


C 

C 


UPDATE THE INDICATOR VECTOR 
vji i = MBASIS ( 


L 

NL2 

MBASIS ( IR) 
MBASIS (L) 
LI 


= N + IR 
= MBASIS ( 
= NE1 
= NE2 
* LI + 1 


OF BASIC VARIABLES 
IR ) 

L ) 


RETURN 

END 


SUBROUTINE SOLVE ( N ) 


PURPOSE - TO CORRELATE COMPLEMENTARY PROBLEM SOLUTION 


C 


C 


c 


c 


c 


c 


c 


COMMON AM( 200 , 200 
1 NE1, NE2, 


) , Q( 200 ) , LI, B ( 200 , 200 ) 
IR, MBASIS( 300 ) , W( 200 ) , 


, NL1 , 
Z( 200 ) 


NL2 


DO 1 I =» I, N 
W(I) = 0.0 

Z(I) = 0.0 

1 CONTINUE 

I - N + 1 

J “ I 


2 K1 
K2 


MBASIS ( I ) 
MBASIS ( J ) 


IF ( Q( J ) -GE. 0.0 ) GO TO 3 
Q(J) = 0.0 

3 IF ( K2 .EQ. 1 ) GO TO 5 
Z(K1) = Q( J ) 

GO TO 7 


5 W(K1) = Q( J ) 

7 1 *1 + 1 
J - J + 1 

IF ( J .LE. N ) GO TO 2 


RETURN 

END 


C 


C 

c**** 

c 

c 


SUBROUTINE PROP ( MAT , PRSI , TMPI , OF , CEFF , ZENTH , ZENTR , 

PROP - PROPERTY PROGRAM CALCULATING HYDROGEN, 


, A( 200 ) , 


ZDENS) 


c OXYGEN, STEAM AND HOT GAS PROPERTIES 

C 

COMMON /H2PRP/ 

* H2P1( 15) , H2T1 (11) ,H2H1{ 15,11) ,H2S1{ 15 ,11) ,H2D1 (15,11), 

* H2P2( 20) ,H2T2( 11 ) ,H2H2( 20 ,11 ) ,H2S2( 20 ,11) ,H2D2(20, 11) , 

* H2P3( 29 ) , H2T3 (25),H2H3(29,25), H2S3 (29,25) ,H2D3( 29,25), 

* H2P4 (23) , H2T4 ( 25 ) ,H2H4(23,25) , H2S4 ( 23 , 25 ) , H2D4 ( 23 , 25 ) 
COMMON /02PRP/ 

* 02P1( 13 ) , 02T1 ( 16 ) , 02H1( 13 ,16) ,02S1(13 , 16 ) , 02D1 ( 13 , 16 ) , 

* 02P2 (13) ,02T2( 17) ,02H2(13,17) ,02S2( 13,17) ,02D2(13,17) , 

* 02P3 ( 5 ) , 02T3 ( 6 1 ) ,02H3(5,61) , 02S3(5,61), 02D3(5,61) 
COMMON /H20PRP/ 

* H20P1 ( 7 ) , H20T1 ( 13 ) ,H20H1(7,13) ,H20S1(7,13) ,H20D1(7,13) 
C 

COMMON /TABLE/ 

* NH2P ( 4 ) , NH2T( 4 ) , N02P ( 3 ) , N02T ( 3 ) , NH20P( 1 ) , NH20T( 1 ) 

C 

DIMENSION 

* TSH2 ( 11 ) , PSH2( 11) , HLH2( 11 ) ,HVH2( 11 ) , SLH2 ( 11 ) , SVH2 ( 11 ) , 

* DLH2( 11 ) , DVH2 ( 11 ) , 

* TS02 ( 16 ) , PS02 ( 16 ) , HL02 ( 16 ) , HV02 ( 16 ) , SL02( 16 ) , SV02 ( 16 ) , 

* DL02( 16) ,DV02( 16) 

C 


c 

TSH2 

- 

H2 

SATURATION 

TEMPERATURE 

c 

PSH2 

- 

H2 

SATURATION 

PRESSURE 

c 

HLH2 

- 

H2 

SATURATION 

ENTHALPY - LIQUID 

c 

HVH2 

- 

H2 

SATURATION 

ENTHALPY - VAPOR 

c 

SLH2 

- 

H2 

SATURATION 

ENTROPY - LIQUID 

c 

SVH2 

- 

H2 

SATURATION 

ENTROPY - VAPOR 

c 

DLH2 

- 

H2 

SATURATION 

DENSITY - LIQUID 

c 

DVH2 

- 

H2 

SATURATION 

DENSITY - VAPOR 


C 

DATA (TSH2 ( J) , J=1 , 11 ) / 

* 30.0,32.0,34.0,36.0,38.0,40.0,42.0,44.0,46.0,48.0,50.0/ 

C 

DATA (PSH2(J) ,J=1,11)/ 

* 4.170,6.446,9.527,13.561,18.694,25.089,32.915,42.334, 

* 53.514,66.625,81.838/ 

C 

DATA (HLH2( J) , J=1 , 11 )/ 

* -123.995,-120.090,-115.893,-111.380,-106.524,-101.289, 

* -95.636,-89.513,-82.850,-75.556,-67.493/ 

C 

DATA ( HVH2 ( J) , J=1 , 11 ) / 

* 70.977,74.584,77.848,80.729,83.256,85.199,86.614,87.431, 

* 87.546,86.817,85.043/ 

C 

DATA ( SLH2 (J) ,J=1,11)/ 

* 1.506,1.629,1.752,1.876,2.002,2.129,2.259,2.391,2.528, 

* 2.670,2.819/ 

C 

DATA ( SVH2 ( J) , J*1 , 11 )/ 

* 8.005,7.713,7.451,7.214,6.998,6.794,6.601,6.415,6.234, 

* 6.054,5.871/ 

C 

DATA (DLH2( J) , J=1 ,11)/ 

* 4.6500,4.5832,4.5127,4.4378,4.3580,4.2724,4.1801,4.0798, 

* 3.9698,3.8479,3.7108/ 

C 

DATA ( DVH2 (J),J=1,11)/ 



* 0.0272,0.0401,0.0568,0.0779,0.1039,0.1363,0.1757,0.2234, 

* 0 . 2809 , 0 . 3508 , 0 . 4362 / 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 


TS02 - 02 SATURATION 

PS02 - 02 SATURATION 

HL02 - 02 SATURATION 

HV02 - 02 SATURATION 

SL02 - 02 SATURATION 

SV02 - 02 SATURATION 

DL02 - 02 SATURATION 

DV02 - 02 SATURATION 


TEMPERATURE 
PRESSURE 
ENTHALPY - 
ENTHALPY - 
ENTROPY - 
ENTROPY - 
DENSITY - 
DENSITY - 


LIQUID 

VAPOR 

LIQUID 

VAPOR 

LIQUID 

VAPOR 


C 

C 

c 


DATA ( TS02 ( J) , J=1 , 16 ) / 

* 160.0,164.0,168.0,172.0,176.0,180.0,184.0,188.0,192.0, 

* 196.0,200.0,204.0,208.0,212.0,216.0,220.0/ 

DATA ( PS02 ( J) , J“1 , 16 ) / 

* 12.810,16.183,20.200,24.935,30.467 ,36.876,44.243,52.654 , 

* 62.194,72.951,85.013,98.473,113.421,129.952,148.162, 

* 168.146/ 

DATA ( HL02 ( J ) , J=1 , 16 ) / 

* -58.356,-56.730,-55.096,-53.455,-51.804,-50.144,-48.473, 

* -46.790,-45.093,-43.380,-41.650,-39.901,-38.130,-36.334, 

* -34.511,-32.657/ 

DATA ( HV02 ( J) , J=1 , 16 ) / ... 

* 33 . 777 ,34.457,35.110,35.734,36.326, 36.884 ,37.408,37.894, 

* 38.340,38.745,39.105,39.419,39.683,39.894,40.049,40.144/ 

DATA <SLQ2(J),J-1,16)/ 

* 0.698,0.708,0.717,0.727,0.736,0.746,0.755,0.764,0.772, 

* 0.781,0.790,0.798,0.806,0.815,0.823,0.831/ 

DATA (SV02( J) , J=l,16)/ , „ , , 

* 1.273,1.263,1.254,1.245,1.237,1.229,1.221,1.214,1.207, 

* 1.200,1.193,1.187,1.180,1.174,1.168,1.162/ 

DATA (DL02(J) ,J-1,16)/ 

* 71. 630 ,70 .941 , 70 • 243 , 69. 536 ,68 • 818 , 68 . 089 , 67 • 347 , 66. 593 , 

* 65.823,65.037,64.234,63.412,62.567,61.699,60.804,59.880/ 

DATA ( DV02 ( J) , J=1 , 16 )/ 

* 0.246,0.305,0.374,0.455,0.547,0.653,0. 774 ,0.911,1.065, 

* 1.239,1.433,1.650,1.893,2.162,2.461,2.794/ 

51 FORMAT (/3X, 'PROP - REQUESTED PRS > ',F7.2,2X, 

* ' AND TMP > ' , F7 . 2 , 2X , ' FOR H2 IS OUT OF RANGE ) 

52 FORMAT (/3X, 'PROP - REQUESTED PRS > ',F7.2,2X, 

* 'AND TMP > ' ,F7.2,2X, 'FOR 02 IS OUT OF RANGE') 

53 FORMAT ( /3X, 'PROP - REQUESTED PRS > '.F7.2.2X, 

* 'AND TMP > ' ,F7.2,2X, 'FOR STEAM IS OUT OF RANGE ) 

** INTERPOLATE RESULTS FROM SINGLE ARRAY ** 

IPRP=0 
NPX1=2 
NPY1=2 
ZENTH=0 . 0 



10 


ZENTR=0 . 0 
ZDENS=0 . 0 

GO TO (10,20,30,40) MAT 

if(tmpi.gt. 30.o. ahd.tmpi.lt. so.o) 

IF(TMPI.GT. 70.0. AND. TMPI.lt. 1 • | IpR p«3 

TrrTMPl GT. 2 4 0 . 0 . AND . TMPI . LT . /20.0) IPRP 

IF TMPl'.GT. 1400.0. AND.TMPI.LT. 2000.0) IPRP=4 

GO TO (11,12,13,14) IPRP 
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O.OR.PRSI.GT. 370.0) GO TO 50 
, TMPI , ZENTH , 

1) , NH2T( 1 ) ,11,29.95,50.05, 

, PSH2 ,TSH2 ,HLH2 ,HVH2 ) 

, TMPI , ZENTR, 

1) , NH2T ( 1 ) ,11,29.95,50.05, 

, PSH2 , TSH2 , SLH2 , SVH2 ) 

,TMPI , ZDENS , 

1) ,NH2T( 1 ) ,11,29.95,50.05, 

, PSH2 , TSH2 , DLH2 , DVH2 ) 


IF(PRSI.LT. 20. 

CALL PRPSAT( PRSI 

* TSH2 (11), NH2P( 

* H2P1 , H2T1 , H2H1 
CALL PRPSAT( PRSI 

* TSH2 ( 11 ) , NH2P( 

* H2P1 ,H2T1 ,H2S1 
CALL PRPSAT( PRSI 

* TSH2( 11) ,NH2P( 

* H2P1 , H2T1 , H2D1 
RETURN 

12 IFfPRSI.LT. 3400. O.OR.PRSI.GT. 7200.0) GO TO 50 
r , T r ITERP2 { PRSI , TMPI , H2P2 , H2T2 , H2H2 , 

* NH2P( 2 ) ,NH2T(2) ,NPX1 ,NPY1 ,NH2P( 2) , ZENTH, Nl) 

CALL ITERP2(PRSI,TMPI,H2P2,H2T2,H2S2,^ > ^ 

* NH2P( 2 ) ,NH2T( 2 ) ,NPX1 ,NFY1 ,NH2P(2), ZENTR, Nl) 

r »Tt TTERP2 ( PRSI , TKPI , H2P2 , H2T2 , H2D2 , 

* NH2P( 2 ) , HH2T( 2 ) ,NPX1 ,NPY1 ,NH2P( 2) , ZDENS, Nl) 

RETURN 

13 IF(PRSI.LT*1400. O.OR.PRSI.GT. 7000.0) GO TO 50 
CALL ITERP2 (PRSI, TMPI, H2P3 ,H2T3 ,H2H3,^^^ 

* NH2P ( 3 ) , NH2T ( 3 ) , NPX1 , NPY1 , NH2P( 3 ) , ZENTH , Nl ) 
CALL ITERP2 (PRSI, TMPI, H2P3,H2T3.H2S3,^ 

* NH2P(3) ,NH2T(3) , NPX1 ,NPY1 ,NH2P( 3 ) , ZENTR, Nl) 

CALL ITERP2 (PRSI, TMPI, H2P3,H2T3,H2D3, 

* NH2P( 3 ) ,NH2T( 3 ) ,NPX1 ,NPY1 , NH2P( 3 ) , ZDENS, Nl) 

return 
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c 

c 


20 


21 


IFfPRSI.LT. 1400. O.OR.PRSI.GT. 5800.0) GO TO 50 
CALL ITERP2 (PRSI, TMPI, H2P4, H2T4 ,H2H4,^^ 

* NH2P(4) ,NH2T(4) ,NPX1 .NPY1 ,NH2P( 4) , ZENTH, Nl) 

CALL ITERP2( PRSI, TMPI, H2P4,H2T4,H2S4,^^ 

* NH2P( 4 ) ,NH2T(4) ,NPX1 ,NPY1 ,NH2P( 4 ) , ZENTR, Nl) 
CALL ITERP2 (PRSI, TMPI, H2P4 ,H2T4 ,H2D4 , 

* NH2P ( 4 ) , NH2T ( 4 ) , NPX1 , NPY1 , NH2P ( 4 ) , ZDENS , Nl ) 

RETURN 

IF(TMPI.GT. 160.0. AND.TMPI.LT. 240.0) IPRP-1 
IF( IPRP.EQ. 1 .AND. PRSI .LT .650 . 0 ) IPRP“1 
IF( IPRP. EQ. 1 .AND. PRSI .GT. 650 . 0) IPRP“2 
IF(TMPI.GT. 600.0. AND.TMPI.LT. 1500.0) IPRP* 

GO TO (21,22,23) IPRP 

IFfPRSI.LT. 30. O.OR.PRSI.GT. 630.0) GO TO 50 
IF TMPI- LT. 160.0. OR. TMPI. GT. 219.9) GO TO 50 
CALL PRPSATf PRSI, TMPI, ZENTH, 



* TS02 ( 16 ) ,N02P( 1 ) , N02T ( 1) ,16,159.95,220.05, 

* 02P1 ,C2T1 , 02H1 , PS02 ,TS02 ,HLC2 , HV02 ) 

CALL PRPSATf PRSI , TMPI , ZENTR , 

* TS02( 16 ) , N02P ( 1 ) , N02T (1) ,16,159.95,220.05, 

* 02P1 , C2T1 ,02S1,PS02,TS02, SL02 , SV02 ) 

CALL ?RPSAT( PRSI, TMPI, ZDENS, 

* TS02 ( 16 > , N02P ( 1 ) , N02T ( 1) ,16,159.95,220.05, 

* 02P1 , 02T1 , 02D1 , PS02 ,TS02,DL02 ,DV02) 

RETURN 

C 

22 IF(PRSI .LT. 2000 . 0 .OR. PRSI .GT. 8000 . 0 ) GO TO 50 
CALL ITERP2 ( PRSI , TMPI , 02P2 , 02T2 , 02H2 , 

* N02P ( 2 ) , N02T( 2 ) , NPX1 , NPY1 , N02P( 2 ) ,ZENTH,N1) 

CALL ITERP2 ( PRSI , TMPI , 02P2 , 02T2 , 02S2 , 

* N02P( 2 ) , N02T( 2 ) , NPX1 , NPY1 , N02P ( 2 ) , ZENTR , N1 ) 

CALL ITERP2 ( PRSI , TMPI , 02P2 , 02T2 , 02D2 , 

* N02P( 2 ) , N02T( 2 ) , NPX1 , NPY1 , N02P ( 2 ) , ZDENS , HI ) 

RETURN 

C 

23 IF(PRSI .LT. 2000. 0 .OR. PRSI .GT. 4000 . 0 ) GO TO 50 
CALL ITERP2 ( PRSI , TMPI , 02P3 , 02T3 , 02H3 , 

* N02P ( 3 ) , N02T( 3 ) , NPX1 , NPY1 , N02P ( 3 ) ,ZENTH,N1) 

CALL ITERP2( PRSI, TMPI, 02P3 , 02T3 , 02S3 , 

* N02P(3) ,N02T(3) , NPX1 , NPY1 , N02P( 3 ) , ZENTR, Nl) 

CALL ITERP2( PRSI , TMPI ,02P3 ,02T3 ,02D3 , 

* N02P( 3 ) , N02T( 3 ) , NPX1 , NPY1 , N02P( 3 ) , ZDENS, Nl) 

RETURN 

C 

C 

30 IF(TMPI.LT. 1400.0. OR. TMPI. GT. 2000.0) GO TO 50 
IF (PRSI .LT. 100.0. OR. PRSI. GT. 700.0) GO TO 50 
CALL ITERP2 ( PRSI , TMPI , H20P1 , H20T1 , H20H1 , 

* NH20P( 1) , NH20T ( 1 ) , NPX1 , NPY1 , NH20P ( 1 ) , ZENTH , Nl ) 

CALL ITERP2 ( PRSI , TMPI , H20P1 , H20T1 , H20S1 , 

* NH20P(1) , NH20T( 1 ) ,NPX1 ,NPY1 ,NH20P(1) , ZENTR, Nl) 

CALL ITERP2 ( PRSI , TMPI , H20P1 , H20T1 , H20D1 , 

* NH20P ( 1 ) , NH20T ( 1 ) , NPX1 , NPY1 , NH20P ( 1 ) , ZDENS , Nl ) 

RETURN 

C 

C 

40 CALL PRPMIX( PRSI, TMPI, OF, CEFF,HMIX,SMIX) 

ZENTH=HMIX 

ZENTR-SMIX 

ZDENS-0.0 

RETURN 

C 

C 

50 IF(MAT.EQ.l) WRITE(21,51) PRSI, TMPI 
IF(MAT. EQ. 2 ) WRITE(21,52) PRSI, TMPI 
IF(MAT. EQ. 3 ) WRITE(21,53) PRSI, TMPI 
RETURN 
C 

END 

C*** ************************* ******************************** 

SUBROUTINE PRPMIX ( P ,TMPI ,OF , CEFF , HMIX , SMIX ) 

C 

C PRPMIX - CALCULATES HOT GAS MIXTURE PROPERTIES. 

C 

COMMON /H2PRP/ 

* H2P1( 15) ,H2T1( 11) , H2H1 (15,11) , H2S1 ( 15 , 11 ) , H2D1 ( 15 , 11 ) , 




* H2P2(20) , H2T2 ( 11 ) ,H2H2( 20,111 , H 2 S2 ( 20 , 11 ) , H2D2 ( 20 , ) , 

* H2P3( 29 ) , H2T3 ( 25 ) ,H2H3( 29 ,25) ,H2S3( 29 ,25) >H2D3(2 . ) * 

* H2P4 ( 23 ) , H2T4 ( 25) ,H2H4 ( 23 , 25) ,H2S4( 23 ,25) , H2D4 (23,2 ) 

* 02P1( 13) ,02T1(16) ,02H1(13 ,16) ,02S1(13 ,16) , 02D1 ( 13 , 16 ) , 

* 02P2 ( 13 ) ,02T2( 17) ,02H2( 13 , 17) ,02S2( 13 ,17) ,02D2(13 ,17) , 

* 02P3 ( 5 ) , 02T3 ( 61 ) ,02H3 ( 5 , 61 ) , 02S3(5,61), 0203(5,61) 
COMMON /H20PRP/ 

* H20P1 ( 7 ) , H20T1 (13) ,H20H1(7,13) ,H20S1(7,13) ,H20D1(7,13) 
COMMON /TABLE/ 

* NH2P ( 4 ) , NH2T( 4 ) ,N02P( 3 ) , N02T( 3 ) ,NH20P( 1) ,NH20T( 1) 
COMMON /STD/ 

* HH2REF,H02REF,HWAREF,SH2REF,S02REF ,SWAREF,SH2A,S02A, 

* SWAA 


XMWH2 = 2.0160 

XMW02 = 31.9988 

XMWH20 = 18.0153 

HCOMB = -6825.6550 

NPX1 = 2 
NPY1 = 2 
ITST1* 0 
ITST2= 0 
ITST3= 0 
ITST4= 0 
ITST5- 0 
ITST6- 0 

XF = 1.0 / (1.0 + OF) 

XO = 1.0 - XF 

XH2 = XF - XO * 2.0 * CEFF * XMWH2 / XMW02 
XH20 = XO * 2.0 * CEFF * XMWH20 / XMW02 

X02 = 1.0 - XH2 - XH20 

EH2 = XH2 / XMWH2 
EH20 = XH20 / XMWH20 
E02 = X02 / XMW02 

ET =■ EH2 + EH20 + E02 

YH2 = EH 2 / ET 

YH20 = EH20 / ET 
Y02 = 1.0 - YH2 - YH20 

PH2 = P * YH2 

PH20 = P * YH20 

P02 = P * Y02 

IF(TMPI .LT. 1000. 0 .OR.TMPI .GT. 2000.0 ) ITST1**1 
IF(PH2 .LT. 1400.0. OR. PH2 .GT. 5800.0) ITST2*1 
CALL ITERP2 ( PH2 ,TMPI , H2P4 , H2T4 , H2H4 , 

* NH2P( 4 ) , NH2T ( 4 ) , NPX1 , NPY1 , NH2P ( 4 ) , HH2 , HI ) 
CALL ITERP2 ( PH2 , TMPI , H2P4 , H2T4 , H2S4 , 

* NH2P ( 4 ) , NH2T( 4 ) , NPX1 , NPY1 , NH2P( 4 ) , SH2 , N1 ) 

IF ( TMPI . LT . 1 4 00 . 0 . OR . TMPI . GT . 2 000 . 0 ) ITST3=1 
IF(PH20.LT. 100 . 0 . OR. PH20.GT. 700.0) ITST4“1 
CALL ITERP2 ( PH20 , TMPI , H20P1 , H20T1 , H20H1 , 
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. NH20P ( 1 ) , NH20T ( 1 ) , NPX1 , ^PYl , NH20P( 1 ) * HH20 , HI 

CALL ITERP2 ( PH20 ,TMPI ,H20P1 ,H20T1 ,H20S < 

‘ NH20P( 1) ,HH20T( 1 ) ,NPX1 ,NPY1 ,NH20P( 1 ) -SH , 


) 

) 


C 


IF( Y02.LT. 0.001) THEN 
DH02 =0.0 
DS02 =0.0 

imMPI.GT. 600.0. AND.THPI.LT. 1500.0) ITST5-1 
I F ( PO 2 . LT. 2000 . 0 .OR. P02 .GT. 4000.0) ITST 

CALL ITERP2(P02,TMPI,02P3,02T3,02H3, 

* H02P ( 3 ) , N02T ( 3 ) ,NPX1 ,NPY1 ,N02P( 3 ) ,H02 ,N1) 
CALL ITERP2( P02 ,TMPI ,02P3 ,02T3 ,02S3 , 

* N02P( 3 ) ,N02T( 3 ) ,NPX1 ,NPY1 ,N02P( 3 ) ,S02,N1) 
DH02 = H02 - H02REF 

DS02 = S02 — S02REF + S02A 
END IF 


C 
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C 


DHH2 * HH2 - HH2REF 
DHH20M =* (HH20 - HWAREF) + HCOMB 
DSH2 = SH2 - SH2REF + SH2A 
DSK20 = SH20 - SWAREF + SWAA 


HMIX 

SMIX 


XH2*DHH2 + XH20*DHH20M + X02*DH02 
XH2*DSH2 + XH20*DSH20 + X02*DS02 


C 

C 


IF ( ITST1 . EQ . I . OR. ITST2.EQ.1) 
IF ( ITST3 .EQ. 1 • OR. ITST4 . EQ. 1 ) 
IF ( ITST5 . EQ . 1 • OR . ITST6 . EQ . 1 ) 


WRITE (21, 51) 
WRITE (21, 52) 
WRITE( 21,53) 


PH2 , TMPI 
PH20,TMPI 
P02 , TMPI 


51 FORMAT ( / 3X , ' PRPMIX - REQUESTED PH2 ep P ^f_ > ' ' pliicEn ' 
* 'AND TMP > ' , F7 . 2 , 2X , ' FOR " H2" *S °UT OF S^^X 

52 FORMAT(/3X, 'PRPMIX - REQUESTED PH20 PRS > ,E7.2,2X, 

“ . ^d™P> ' , F7 . 2 ,2X, ' FOB "H20" IS OUT OF RAHOE') 

„ FQRMATf /3X. 'PRPMIX - REQUESTED P02 PRS > ,F7.2,2X, 

” / ^ > ^F7.2,2X,'FOR " 02" IS OUT OF RANGE') 


return 

END 


SUBROUTINE PRPSAT ( X , Y , FPROP , TCRT , NX1 , NY1 , NX2 , YL , YH , 

* PRS1 , TMPI , PROP , PRS2 , TMP2 , PROPL , PROPV) 

PRPSAT “ CALCULATES NBS PROPERTIES NEAR SATURATION CURVE 


DIMENSION PRS1 ( 1 ) , TMPI ( 1 ) 
C 

NR1-NX1 

NPX1~2 

NPY1=2 

NPX2=2 

C 

ZPLGAS-0 . 0 
ZPHGAS-0.0 
ZPLLIQ-0.0 
ZPHLIQ=0 . 0 
ZPR0P1“0.0 
ZPROP-O . 0 
FPROP=0 . 0 
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ZTSAT=0 . 0 
ARGA=0 . 0 
AJRGB=0 . 0 
ZTSAT7=0 . 0 


CALL ITERP2 ( X , Y , PRS1 , TMP1 , PROP , NX1 , NY1 , NPX1 , NPY1 , NR1 , ZPROP1 , N1 ) 

FPROP=ZPROP 1 

IF ( Y.GT.TCRT) GO TO 70 

CALL ITERP1 ( X , PRS2 , TMP2 , NX2 , NPX2 f ZTSAT , N2 ) 

IF( Y. LT. ZTSAT) GO TO 61 

* * GAS CALCULATIONS * * 

CALL ITERP1 ( X , PRS2 , PROPV , NX2 , NPX2 , ZPGAS , N2 ) 

CALL ITERP2 ( X , YH , PRS 1 , THP1 , PROP , NX1 , NY1 , NPX1 , NPY1 , NR1 , ZTST f N1 ) 

DTST- ZTST- ZPGAS 

IF(DTST.GT. 0.0001) GO TO 50 

ZPLGAS= ZPGAS 

IF (ZPROP1.LT. ZPGAS) GO TO 70 
GO TO 51 

50 ZPHGAS=ZPGAS 
IF(ZPR0P1.GT. ZPGAS) GO TO 70 

51 LPR*1 

53 PRSD=PRS1(LPR) -0.0001 
IF(PRSD.GT.X) GO TO 52 
LPR-LPR+1 

GO TO 53 

52 ARGA=PRS1 ( LPR) 

CALL ITERP 1 ( ARGA , PRS 2 , TMP2 , NX2 , NPX2 , ZTSATT , N2 ) 

LTP=1 

54 TMPD=TMP1(LTP) -0.0001 
IF(TMPD.GT. ZTSATT) GO TO 55 
LTP-LTP+1 

GO TO 54 

55 ARGB-TMPl(LTP) 

YY*ARGB 

IF(DTST.GT. 0.0001) CALL ITERP2 ( X , YY , PRS1 , TMP1 ,PROP , NX1 # NY1 , 

* NPX1 ,NPY1 , NR1 , ZPLGAS r N1 ) 

IFfDTST.LT. 0.0001) CALL ITERP2(X,YY,PRSl f TMPl f PROP,KXl ,NY1 , 

* NPX1 f NPY1 , NR1 , ZPHGAS , N1 ) 

ZPROP=ZPHGAS-( ZPHGAS- ZPLGAS) * ( (ARGB-Y)/(ARGB- ZTSAT) ) 

FPROP*ZPROP 

GO TO 70 

* * LIQ CALCULATIONS * * 

61 CALL ITERP1 ( X , PRS2 , PROPL , NX2 f NPX2 , ZPLIQ , N2 ) 

CALL ITERP 2 ( X, YL, PRS1 , TMP1 , PROP , NX1 , NY1 , NPX1 , NPY1 , NR1 , ZTST , N1 ) 

DTST= ZTST- ZPLIQ 

IF(DTST.GT. 0.0001) GO TO 59 

ZPLLIQ= ZPLIQ 

IF(ZPR0P1.LT. ZPLIQ) GO TO 70 
GO TO 60 
59 ZPHLIQ=ZPLIQ 


IFCPRCP1.GT.ZPLIQ) go TO 70 

c 

60 LPR = 1 

63 PRSD=FRS1( LPR ) -0 . 000 1 
IF(FRSO.GT.X) GO TO 62 
LPR=LFR*1 

GO TO 63 
C 

62 ARGA=?RS1 ( LPR-1 ) 

CALL ITERP1 ( ARGA , PRS2 , TMP2 , NX2 , NPX2 , ZTSATT , N2 ) 

C 

LTP=1 

64 TMPD=TMP1( LTP) -0.0001 

IF (TMPD.GT. ZTSATT) GO TO 65 
LTP=LTP+1 
GO TO 64 
C 

65 ARGB=TMP1 (LTP-1 ) 

YY=ARGB 

IF(DTST.GT. 0.0001) CALL ITERP2 ( X , YY , PRS1 ,TMP1 r PROP ,NX1 , NY1 , 

* NPX1,NPY1,NR1,ZPLLIQ,N1) 

IF(DTST.LT. 0.0001) CALL ITERP2 ( X , YY , PRS1 ,TMP1 ,PROP , NX1 ,NY1 , 

* NPX1 , NPY1 , HR1 , ZPHLIQ , N1 ) 

ZPROP=ZPHLIQ- ( ZPHLIQ-ZPLLIQ ) * ( ( ZTSAT-Y ) / ( ZTSAT-ARGB ) ) 
FPROP=ZPROP 
C 

70 CONTINUE 
C 

RETURN 

END 

c******* ********* *★*★*★****★***★ ********************************** 

SUBROUTINE ITERP1 ( X , XT , YT , NX , NPX , Y , NERR) 

C 

C ITERP1 - SINGLE INTERPOLATION ROUTINE. 

C 

DIMENSION XT ( 1 ) , YT ( 1 ) 

NERR=Q 

INTER=1 

NP-NPX 

IF (NX .LT. NP) NP=NX 

IH=NP/2 

1=1 

IF(XT( I )-X) 30,20,10 
10 IH-0 

12 NERR=1 
GO TO 70 

13 NERR=2 
GO TO 70 

20 INTER=2 
22 Y«YT(I) 

GO TO 999 
30 I=NX 

IF(XT(I)-X)13 ,20,40 
40 Nl-1 
N2=NX 

45 MP=(Nl+N2)/2 

50 IF(XT(MP)-X)52,54,56 

52 N1=MP 

GO TO 60 
I=MP 


54 
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GO TO 20 
56 N2=KP 

60 IF! ( N2-N1 ) .NE. I) GO TO 45 

IF (N2.GT. (IH+1) ) GO TO 65 
I=IH-1 
GO TO 70 
65 I=N2 


IF ( N 2 .GT. I) I=N2 
70 K=I-IH 

N=K+HP-1 
Y=0 . 

IF(N-NX)90,90,80 
80 N=HX 

K=NX-NP+1 
90 DO 120 J=K , H 
P=1.0 

DO 110 I=K, H 
IF(I-J)100, 110,100 
100 P=P* ( X-XT( I ) )/(XT( J)“XT( I) ) 
110 CONTINUE 

Y=Y+YT( J)*P 
120 CONTINUE 
GO TO 999 

ENTRY ENTERP (X,XT,YT,Y) 
Y»0. 

GO TO (90,22) , INTER 
999 CONTINUE 
RETURN 


END 

C******** 


**************************************************nnt** 


SUBROUTINE ITERP2 (X,Y,XT,YT,ZT,NX,NY,NPX,NPY,NR,Z,NERR) 


ITERP2 - DOUBLE INTERPOLATION ROUTINE. 


DIMENSION XT( 1 ) , YT( 1 ) , ZT(NR, 1 ) , ZC( 15 ) 

NERRB=0 

NPYY=NPY 

IF (NY .LT. NPY) NPYY=*NY 
IH=NPYY/2 
1=1 

IF(YT(I)-Y)30,20,10 
10 IH=0 

12 NERRB-201 
GO TO 70 

13 NERRB=204 

GO TO 70 , . 

20 CALL ITERP1( X, XT, ZT( 1 , 1 ) ,NX ,NPX, Z ,NERRA) 
GO TO 999 
30 I=NY 

IF( YT( I ) -Y ) 13 , 20 , 40 
40 Nl=l 
N2=NY 

45 MP=( N1+N2 ) /2 
50 IF(YT(MP)-Y)52,54 ,56 
52 N1=MP 
GO TO 60 
54 I=MP 

GO TO 20 



56 N2=MP 

60 IF( ( N2-N1 ) .NE. 1) GO TO 45 
I=N2 

IF( I .LT. ( IH+1 ) ) I=IH+1 
70 K=I-IH 

N=K+NPYY-1 
IF(N-NY)90,90,80 
80 N=NY 

K=NY-NPYY+1 
90 J=0 

DO 100 I=K, N 
J=J+1 

IF( J .NE. 1) GO TO 95 

CALL ITERP1 ( X , XT , ZT ( 1 , 1 ) , NX , NPX , ZC ( J ) , NERRA) 
GO TO 100 

95 CALL ENTERP ( X , XT , ZT ( 1 , 1 ) ,ZC(J) ) 

100 CONTINUE 

CALL ITERP1 ( Y , YT(K) , ZC , NPYY , NPYY , Z , NERRC ) 

999 NERR-NERRA+NERRB 
RETURN 
END 
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